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ABSTRACT 

The Cassini flyby of Jupiter in 2000 provided spatially resolved spectra of Jupiter’s atmosphere 
using the Visual and Infrared Mapping Spectrometer (VIMS). A prominent characteristic of these 
spectra is the presence of a strong absorption at wavelengths from about 2.9 jim to 3.1 /im, previously 
noticed in a 3-/im spectrum obtained by the Infrared Space Observatory (ISO) in 1996. While Brooke 
et al. (1998, Icarus 136, 1-13) were able to fit the ISO spectrum very well using ammonia ice as the sole 
source of particulate absorption, Sromovsky and Fry (2010, Icarus 210, 211-229), using significantly 
revised NH 3 gas absorption models, showed that ammonium hydrosulfide (NH 4 SH) provided a better 
fit to the ISO spectrum than NH 3 , but that the best fit was obtained when both NH 3 and NH 4 SH were 
present in the clouds. Although the large FOV of the ISO instrument precluded identification of the 
spatial distribution of these two components, the VIMS spectra at low and intermediate phase angles 
show that 3-/tm absorption is present in zones and belts, in every region investigated, and both low- 
and high-opacity samples are best fit with a combination of NH 4 SH and NH 3 particles at all locations. 

The best fits are obtained with a layer of small ammonia-coated particles (r ~ 0.3 fim) overlying but 
often close to an optically thicker but still modest layer of much larger NH 4 SH particles (r ~ 10 
/rm), with a deeper optically thicker layer, which might also be composed of NH 4 SH. Although these 
fits put NH 3 ice at pressures less than 500 mb, this is not inconsistent with the lack of prominent 
NH 3 features in Jupiter’s longwave spectrum because the reflectivity of the core particles strongly 
suppresses the NH 3 absorption features, at both near-IR and thermal wavelengths. Unlike Jupiter, 

Saturn lacks the broad 3-/rm absorption feature, but does exhibit a small absorption near 2.965 /tm, 
which resembles a similar Jovian feature and suggests that both planets contain upper tropospheric 
clouds of sub-micron particles containing ammonia as a minor fraction. 

Subject headings: Jupiter; Jupiter, Atmosphere; Jupiter, Clouds 


1. INTRODUCTION 

Analysis of Pioneer a nd groundbased o bservations of 
Jupiter, summarized bv lWest et al.l (jl986D . led to an ex¬ 
pected Jovian cloud structure that included an upper 
ammonia cloud layer starting near 700-mb and a puta¬ 
tive NH 4 SH cloud top near 2 bars, which was thought to 
be optically thick outside the hot spot regions. The pu¬ 
tative ammonia cloud was thought to have two particle 
populations: a vertically compact layer of large particles 
(of 3 to 100 fim in radius) and a vertically diffuse com¬ 
ponent of small particles (r ~ 1 /rm) extending up to 200 - 
300 mb in low latitude regions. However, the more diffuse 
component should have produced prominent spectral sig- 
natures at 9.4 /rm and 26 /rm, which were not observed 
(|Orton et al.lll98^ . After considering possible masking 
of these f eatures by the likel y tetrahedral shapes of these 
particles, IWest et HI (|1989l) concluded that these parti¬ 
cles could not be primarily composed of ammonia ice. It 
was also the case that neither Voyager Inf rared Interfer¬ 
omet er Spectrometer (IRIS) observ ations (iCarlson et all 
II 993 H , nor microwave observations (|de Patedll986[) . ever 
found an ammonia vapor profile that would support a 
700-mb condensation level. Carlson’s derived condensa¬ 
tion pressure was closer to 500 mb, while a later analy- 

^ University of Wisconsin - Madison, Madison WI 53706 
t Partly based on observations obtained from the data archive 
at the Space Telescope Science Institute. STScI is operated by 
the Association of Universities for Research in Astronomy, Inc. 
under NASA contract NAS 5-26555. 


sis of microwave observations (|de Pater et al.l[20MI) sug¬ 
gested NH 3 condensation near 600 mb. 

Only in the last decade or so has there been even a hint 
of the spectral signatures expected of NH 3 ice clouds. 
From an analysis of a 3-/im a bsorption anomaly i n a 
central-disk spectrum of Jupiter. [Brooke et al.l (119981) in¬ 
ferred the existence of a layer of ammonia ice particles 
of 10 fim in radius, beginning at 550 mb with a scale 
height of 30% of the gas scale height. The wide field of 
view covered by the observation (roughly a quarter of the 
Jovian disk) suggested that the ammonia ice was widely 
distributed, lirwin et all (|2001D found a similar absorp¬ 
tion anomaly in analysis of observations by the Galileo 
Near Infrared Mapping Spectrometer (NIMS), but con¬ 
cluded that it was not due to NH 3 ice because a key spec- 
tral signatur e at 2. 0 fiia was missing. The subsequent 
I Baines et al.l (I2002D detection of spectrally identifiable 
ammonia clouds (SIACs) in other NIMS observations was 
based on depressed reflectivity at 2.7 fim as well as at 
2.0 /im, but these detections covered a very tiny fraction 
(< 1%) of Jupiter ’s cloud features. On the other hand, 
iWong et al.l (120041) inferred more widely distributed am¬ 
monia ice, at least in some latitude bands, from a detec¬ 
tion of a 9.4-/im spectral feature in Jovian spectra ob¬ 
tained from the Cassini Composite I nfrared Spec t rome- 
ter (CIRS). The model calculations of IWong et al.l (I2004D 
implied that the 9.4-/rm ice feature could only be de¬ 
tected when the aerosols were at pressures <50(1 mb and 
the particle effective radius was within a factor of two of 
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New band mode ls for NH 3 absorption developed by 
iBowles et al.l (|2008f) and r ecent additions to the HITR AN 
line data base motivated iSromovskv and Fr^ (j20ir)fl to 
update exponential sum models for NH 3 absorption. 
They found substantial changes rel ative to the absorp - 
tion models used in the analysis bv iBrooke et al.l (119981 1 
and undertook a reanalysis of the ISO spectrum. They 
found that the best fit to the 3-^m absorption feature was 
obtained not with NH 3 as the sole absorber, but instead 
when particles of NH 4 SH and NH 3 were both present, 
with the predominant opacity provided by NH 4 SH. Two 
types of solutions were found to provide good fits to the 
ISO spectrum: one had relatively large NH 3 and NH 4 SH 
particles at nearly the same pressure, near 500 mb, and 
the other had small NHs-coated particles closer to 350 
mb, with large NH 4 SH particles again near 500 mb. The 
large FOV of the ISO measurement did not permit any 
conclusions concerning the spatial distribution of these 
two components, i.e. whether some clouds had only 
NH 4 SH particles, and others had only NH 3 particles, or 
whether the two types existed in the same physical loca¬ 
tion. 

In the following analysis we present new results based 
on near-IR spatially resolved spectral observations of 
the Cassini Visual and Infrared Mapping Spectrometer 
(VIMS). Our analysis takes advantage of recently im¬ 
proved models of methane and ammonia gas absorption, 
and uses refractive index properties of putative cloud ma¬ 
terials in constraining model structures and composition. 
We first describe the VIMS observations, then the radi¬ 
ation transfer methods, the modeling results, how they 
compare with other models based on other observations 
at visible, near-IR, and thermal wavelengths, and finally 
summarize their implications for Jovian cloud structure. 
We show that VIMS observations display a widely dis¬ 
tributed strong absorption in the 3.1-/xm region of the 
spectrum, and evidence for a much weaker absorption 
near 2 /xm, both of which are qualitative features ex¬ 
pected from a cloud consisting of ammonia ice particles. 
However, it is only the darker cloud features, with smaller 
optical depths, that can be well fit with ammonia ice as 
the sole absorber. We find that brighter clouds require 
more absorption at 3 /xm without increased absorption 
at 2 /xm and that this can be provided by NH 4 SH. 

2. OBSERVATIONS. 

2 . 1 . VIMS Overview 

For constraining models of low-latitude Jovian cloud 
structure we used spectra acquired by VIMS during 
Cassini’s flyby of Jupiter in December 2000. The 
VIMS instrument and investigation are described by 
iMiller et al.l (|1996f l and iBrown et al.l (j2004[ l . The instru- 
ment’s spectral range is 0.35-5.1 /xm, with an effective 
pixel size of 0.5 milliradians on a side and a near-IR 
spectral resolution of approximately 15 nm (sampled at 
intervals of approximately 16 nm). Our analysis focuses 
on the 1-3.2 /xm interval where scattering by aerosols 
predominates over Rayleigh scattering and where a wide 
range of atmospheric pressures can be sampled and dis¬ 
criminated by means of variations in gas absorption with 
wavelength (Fig.[T]). In this wavelength range the contri¬ 
bution of Jupiter’s thermal emission is everywhere less 


than 0.1% of the observed I/F. To provide an approx¬ 
imate comparison with prior low phase angle observa¬ 
tions we chose VIMS data cube CM_1355256529, ob¬ 
tained on 11 December 2000 at a phase angle of 2.5° 

. This provides a pixel size of 9925 km at the sub¬ 
spacecraft point. In Fig. [2] VIMS images extracted from 
this cube at 4.798 /xm, 1.864 /xm and 1.997 /xm are com¬ 
pared to a Cassini ISS (Imaging Science Subsystem) im¬ 
age obtained at 0.451 /xm. The encircled pixels in the fig¬ 
ure are locations from which we extracted VIMS spectra, 
taking samples from the North Equatorial Belt (NEB), 
and from the Equatorial Zone (EZ). The pixel labeled 
12xl6y, which is brightest at 4.798 /xm, is located in a 
region where deep cloud opacity is at a relative minimum, 
corresponding to a relatively low reflectivity at 1.864 /xm. 
To obtain samples at higher spatial resolution and less 
sensitive to details of backscatter phase functions, we 
also selected data cube CM_1356976257, obtained on 31 
December 2000 at a phase angle of 67.8°, with a sub¬ 
spacecraft pixel size of 4887 km, which is about twice 
the resolution obtained from the low phase angle cube. 
Sample images and target selections from this cube are 
shown in Fig. [H 

2 . 2 . Image Processing and Navigation 

VIMS image cubes were processed using the pipeline 
processing code dow nloaded from the Pla netary Data 
System, described bv IMcCord et al.l (|2004ll . The cubes 
were navigated by manual adjustments of planet center 
coordinates using visual comparisons of model and im¬ 
age limbs at window wavelengths to find an approximate 
best fit, which is estimated to be within a fraction of a 
pixel. This provided sufficient definition of view angles 
and latitudes near the disk center to allow a meaningful 
comparison of spectral characteristics in broad latitude 
regions, as appropriate for the relatively low spatial res¬ 
olution of these observations. 

2.3. Photometry 

The VIMS calibration obtained by the pipeline pro¬ 
cessing has an uncertainty that has not yet been doc¬ 
umented, but is at least as large as the near-IR solar 
spect ral irradiance uncertainty of ~3%-5% (|Colina et al.l 
II996I1 . It appears to be less than 10% uncertain, based 
on window-region comparisons with independent spec¬ 
tra of Jupiter, shown in Fig. 21 In the 0.9-0.98 /xm re¬ 
gion the disk-i ntegrated VIMS sp ectrum is ^ 8 % brighter 
than that of iKarkoschkal (|1998f l. In the 0.9-2.2 /xm 
region we computed a central-disk VIMS spectrum to 
compare with the ground based central-disk spectrum of 
iClark and McCordl (1197911 . In the brighter window re¬ 
gions, where thin upper-level hazes and low-level arti¬ 
facts have much reduce d impact, the VIMS observ ations 
are within ~5% of the IClark and McCordl (|1979ll spec¬ 
trum, except near 1.6 /xm, where an order-sorting filter 
joint in the VIMS detector array apparently results in 
poor response corrections. The VIMS spectrum aver- 
aged over the same a ngular FOV as a 1996 ISO spectrum 
(|Brooke et al.lll99^ is found to be about 20% brighter, 
which might be entirely due to temporal variations or 
even phase angle variations (the VIMS and ISO phase 
angles in this comparison are 2.5° and 11° respectively). 
Another possible contribution to the VIMS - ISO dif¬ 
ference is the ISO calibration uncertainty, which could 
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Fig. 1.— Penetration depth of photons in a clear Jovian atmosphere, indicated by pressures at which a unit change in albedo of a 
reflecting layer results in an external I/F change of 1/e, 1/e^, and l/e‘^ (three solid curves), for NH 3 mixing ratios typical of the NEB. The 
sloping dotted lines define the location at which 2-way optical depths for Raleigh scattering reach 0.1 and 0.01. 



Fig. 2.— Ceissini VIMS images from 11 December 2000 at 4.798 fim (A), 1.864 ^m (B), and 1.997 ^m (C) compared to an ISS image 
acquired roughly one month earlier (D). Circles in A-C, locate pixels from which we extracted spectra for analysis. These are located at 
6 °N, 290°E (12xl6y), 1.2°S, 306°E (llxl4y), and 14°N, 306° E (13xl4y). An additional spectrum (not shown) was extracted from 1.5°S, 
298°E (llxl5y). Here y increases to the left and x upward. 



Fig. 3.— Cassini VIMS images from 31 December 2000 at 4.798 /i.m (A), 1.864 fim (B), and 1.997 fim (C). Circles in A-C, locate pixels 
from which we extracted spectra for analysis. These are located at 9.9° N (centric), 298° E (25x,35y), 2.4°N, 295°E (27x34y), and 1.8°S, 
280°E (28x31y). An additional spectrum was extracted at 6.3°N, 280°E (26x31y). Here y increases to the right and x downward. 


be as large as 12% for wavelengths less than 3.8 dency for the VIMS spectra to be on the high side, but 

(|Schaeidt et al.l[l996f) . Overall, there seems to be a ten- how much of the differences are due to calibration errors 
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as opposed to temporal variation or phase-angle differ¬ 
ences is uncertain. 

2.4. VIMS artifacts. 

There are two important artifacts that need to be con¬ 
sidered in the interpretation of VIMS spectra. The first 
has to do with responsivity corrections at wavelengths 
in close proximity to joints between order sorting filters 
overlaid on the VIMS linear detector arrays. The joints of 
concern here appear at wavelengths of 1.64 /rm and 2.98 
fim. The required responsivity correction factors to com¬ 
pensate for shading effects are very large for the former 
(~300) and relatively small for the latter (~2). To test 
the correction near 2.98 fim w e used a comparison be¬ 
tween the ISO 3-^m spectrum of iBrooke et al.l (|1998D and 
a central disk VIMS spectrum covering approximately 
the same field of view. When the ISO spectrum is em¬ 
pirically scaled by a factor of 1.2, the two spectra are 
in generally excellent agreement (Fig. [SJ, including the 
region of the filter joint near 2.98 /im. However, we do 
see a significant disagreement near 2.58 /rm, where the 
ISO spectrum reaches a deeper minimum. This may be 
related to the second kind of artifact, which is discussed 
later. There is also a small discrepancy near 2.52 /im, 
where the ISO spectrum has a local minimum that is not 
seen in the VIMS spectrum. A similar discrepancy is also 
seen between model spectra and the VIMS observations 
in this region. The rather large responsivity correction 
near 1.64 /im does have significant errors, which can be 
se en from a comp a rison b etween a groundbased spectrum 
of iBanfield et al.l (|1998all and a VIMS spectrum from a 
nearby latitude region (Fig. IH]). The discrepant region is 
from about 1.60 to 1.68 um. This can also be seen from 
a comparison to the iClark and McCordI (|1979t) spectrum 
shown in Fig. |4p. It is also apparent in the striking dis¬ 
agreement between model and observations in this re¬ 
gion, as discussed in later sections. 

The second type of artifact is evident in the 2.3-/tm 
spectral region where the groundbased spectrum in Fig. 
[6] has a minimum that is a factor of 10 lower than that of 
the VIMS spectrum. This cannot be attributed to local 
variations because the latitude difference is small and the 
I/F disagreement is typical of several comparisons. To 
resolve which of these two measurements is more believ¬ 
able we used 1997 NICMOS observations (HST Program 
7445, Reta Beebe PI), selecting an image (n49o01mvs) 
made with the F237M filter, which heavily weights the re¬ 
gion of disagreement, as shown by the filter transmission 
plot in Fig. [S] We computed synthetic F237M I/F val¬ 
ues from the groundbased I/F spectrum and from VIMS 
I/F spectra along a central meridian scan using weighted 
integrals over the spectra. The weighting function was 
the product of filter throughput and solar spectral irra- 
diance, normalized to yield a unit integral over the same 
spectral range. As sh own in Fig. [71 t he synt hetic pho¬ 
tometry result for the IBanfield et al.l (Il998al) spectrum 
is in good agreement with the NICMOS scan. But syn¬ 
thetic F237M values obtained from VIMS spectra are 
^4.5 times as large as the NICMOS values. The artifact 
is not a simple offset error, because that does not provide 
the correct spectral shape. It seems more likely that the 
artifact is due to light scattered by the grating, which 
allows light in bright regions of the spectrum to contam¬ 
inate measurements of dark regions. This may also be 


responsible for the discrepancy with the ISO spectrum 
near 2.58 /rm evident in Fig. [5l In any case, the defect 
in the VIMS measurements near 2.3 /rm has the effect of 
distorting our model fits in the stratosphere, suggesting 
far more particulate reflectivity than is actually present. 

2.5. VIMS noise characteristics. 

For most applications at moderate signal levels, the 
random noise level of VIMS is very small, nominally 
less than one digital quantization number (DN) (Kevin 
Baines, personal communication, January 2010). As 
there is no published equation defining how VIMS mea¬ 
surement noise depends on signal level, exposure time, 
and wavelength, we tried to better characterize the noise 
by comparing spectra in cubes that imaged the same re¬ 
gion of Jupiter with as small a time difference as possible. 
We chose CM_1355309574 and CM_1355309697 cubes, 
which were exposed only two minutes apart. These re¬ 
sults were consistent within about 1.2 DN RMS for a 
single measurement, even for signal levels up to 2000 DN 
or more. For the cubes we used in our detailed spec¬ 
tral analysis, which had exposures much shorter than 
these comparison cubes, the noise level would still be 
only 0.14% at 2.7 /im and thus not a significant error 
source. However, at signals 100 times smaller, as found 
near 2.3 /im, the noise becomes a much larger fraction 
of the total signal. This is evident from the rapid varia¬ 
tions with wavelength observed in that region of the spec¬ 
trum, which are up to 50% or more of the signal level. 
This can be roughly characterized as an offset noise in 
I/F, which we crudely estimated as ^5x10“"^ from the 
on-disk spectra. Other possible sources of uncertainty, 
though not strictly random noise, are scattered light in¬ 
side the spectrometer, wavelength errors, line-spread un¬ 
certainties, and variable absolute calibration and wave¬ 
length calibration errors. We lumped all these potential 
uncertainties into a somewhat arbitrary 2.5% error pro¬ 
portional to signal, and an I/F offset error of ~5xl0“^. 
We did look into the possibility of wavelength errors by 
making fits with shifted wavelength scales, and determin¬ 
ing as a function of shift. We found that the minimum 
was within 1 nm of the nominal wavelength scale and 
that no revision was needed. 

3. RADIATION TRANSFER CALCULATIONS 
3.1. Atmospheric structure and composition 

We followed th e sam e approach described by 
iSromovskv and FrvI (1201011. In sum mary, we used the 
tabulated results of iSeiff et al.l (I1998D for Jupiter’s tem¬ 
perature structure down to the 22-bar level, and assumed 
an atmospheric comp osition of He/H2=0.157±0.003 
( von Zahn et~aD [19^ and CH 4 /H 2 = 2.1E-3±0.4E-3 
( Niemann et al.l 1199811 . which are expressed as number 
density ratios. We used model mixing ratio profiles (Fig. 
0 ) ranging from what should be typical of the equatorial 
zone to hot- spot values. The basis for these models is 
discussed bv iSromovskv and Fr-d (1201611 . The deep mix¬ 
ing ratio is from dFM cner et al.||199 8fl. the u pper level 
Probe result is from ( Sromovskv et all 119981) . The EZ 
(Equatorial Zone, ±7° latitude), NEB (North Equatorial 
Belt, r; 8°N-18°N), Hot Spot (ft!7°N), and EZG ( ±7° lat¬ 
itude) profiles were based on combinations of NH,^ mixing 
ratios inferred from ISO observations by iFouchet et al.l 
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Fig. 4.— VIMS I/F spectra (solid lines) compared to the IKajj^schl^ ([19^) 1995 geometric albedo spectrum of Jupiter (A), the 
IBrooke et al.l 11199^1 ISO 1996 spectrum (B), and the 1976 IClark and McCordl II1979I) central-disk spectrum (C). In (A) and (B) the 
reference spectra (dot-dash) are averaged to VIMS resolution. In (C) the VIMS spectrum (solid) is averaged to the resolution of the 
reference spectrum (dot-dash). Ratio plots are dotted, except where heavy black lines are used in low-absorption regions that are not much 
affected by upper level aerosols. 



Fig. 5.— Comparison of central-disk I/F spectra from VIMS 
(solid) and ISO x 1.2 (dot-dash) with comparable fields of view cov¬ 
ering about 1/4 of the disk. The ISO spect rum was converted from 
the flux spectrum of IBrooke et al.l 1)19981') . Vertical dotted lines 
mark the locations of an absorption feature of NH 3 ice at 2.95 /rm 
and the VIMS filter joint at 2.98 /rm. 

(1200^ and on mi c rowav e results of iSault et al.l (|2004[1 
and iGibson et all (|2005ll . Upper level mixing ratios 


were also guided by r esults of lEdging ton et al.l (jl999H . 
iLara et alJ ()1998f ). an d lAchterberg et alJ (j2006f l . 


3.2. Gas absorption models 

We used the k-distribution models of hydrogen- 
broad ened methane absorption described bv lirwin et al.l 
(|2006ll . which em ploy new two-term mod els of tempera¬ 
ture dependence (iSromovskv et aT1l2006fl . For ammonia 
we used t he combined correlated-k ab sorption model de¬ 
scribed bv ISromovskv and FrvI (l2ni(Hl. which is based o n 
the Goody-Lorentz band model of iBowles et al.l (|2008ll . 
and also makes use of 1996 and 2004 (or equivalently, 
2008) HITRAN spectral line compilations, he basis for 
and detail s of how these resul t s wer e combined are de¬ 
scribed by [Sromoysk^ani^^ (HMS) • Collision-induced 
absorption (CIA) for H 2 and H 2 -He was calculated us¬ 
ing programs downloaded from the Atmospheres Node 
of the Planetary Data System, which are documented by 
Borvsowl lll99lL II 993 II for t he H 2 -H 2 fundamental band, 
Zheng and Borvsow (1199511 fo r the first H 2 -H 2 overtone 
band, and bv iBorvsowl (Il99ll for H 2 -He bands. We fol- 
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Fig. 6 .— Top: Comparison of VIMS 6 ° N sp ectrum (solid) 
with a 7° N groundbased spectrum (dot-dash) from IBanfield et ahl 
Ill998al) . Vertical dotted lines at 1.6 and 1.68 fim bound the region 
in which VIMS spectra are in significant error. Lower: transmis¬ 
sion spectrum of the NICMOS F237M filter. 



Centric Latitude, ° 

Fig. 7. — Comparison of NICMOS F237M central meridian I/F vs 
latitude (solid) with synthetic photometry band-pa ss results com- 

J rted f rom VIMS spectra and a 7° N spectrum from IBanfield et al.1 
-998al l. The filled circles are from VIMS locations identified in Fig. 
The plus symbols are for a central meridian latitude scan with 
a 3-pixel longitudinal average. Offset (dashed) and scaled (dotted) 
versions of the NICMOS scan are also shown. 

lowed the iBirnbaum et al.l (|1996D parameter recommen¬ 
dations to avoid sudden drops in wing absorption in the 
1.58-/tm window. We assumed equilibrium mixing ratios 
of ortho and para hydrogen in all our model fitting cal¬ 
culations and later show that the resulting errors are not 
significant to our main conclusions. 

3.3. Reflecting layer methods 
3.3.1. Linear combination equations 

If all the cloud layers were composed of broken cloud 
fields of opaque elements, a linear combination of reflect¬ 
ing layer contributions would be an ideal model for repre¬ 
senting the combined I/F of the vertical cloud structure. 
By fitting parameters of such a model to obtain agree¬ 
ment with the observed I/F spectra or variation of I/F 
with angle, it can also be used to obtain a rough idea of 
the distribution of scatterers when the cloud layers are 



Fig. 8. — Model NH 3 profiles, adapted from ISromovskv and Fr^ 
(goTO). All profiles overlap above the point near 330 mb where they 
intersect the Galileo Probe profile (lower solid curve) and below 
the intersection near 4.2 bars. At intermediate pressures the Hot 
Spot (dot-dash), NEB (triple-dot dash), EZ (dash-h), and EZG 
(dotted-h) follow separate paths until intersecting the saturation 
mixing ratio (upper solid curve), at which point all but the EZG 
profile follow the same line to the upper intersection with the Probe 
profile. 


optically thin. Its main virtue is the great speed of the 
inversion, because it makes use of precomputed arrays. 
Our model equat ions for a multilaye r atmosphere are 
based on Eq. 4 of iSromovskv and FrvI (|2007D . which can 
be written in the following form in a three-layer model: 


/(0)=/o(A,0) + /l X [Ji(Pi,0)-/o(Pi,0)] 


+ 


/o(^2,0)-/o(A,0) + /2 X 


[h{P2,e) - h{P2,e)] xTflO) 


+ 


/o(P3,0)-/o(P2,0) + /3 X 


[h{P3,e) - h{P3,e)] xTfl9)T2{9) 


+ 


h{9) - loiPs) X Tfl9)T2{9)Tfl9) (1) 


where /o(P, 9) is the I/F observed at the top of the atmo¬ 
sphere for viewing geometry 9 when a zero-albedo sur¬ 
face is placed at pressure P and Ji(P, 0) is the top-of- 
atmosphere I/F for a unit-albedo surface placed at the 
same level. Here 9 denotes the viewing geometry vector 
[0, 00 , (/], which refer to viewer and solar zenith angles 
and azimuth angle, respectively. 

The first line in the above equation contains one con¬ 
tribution from the atmosphere above the first reflecting 
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layer (/q) and a second differential contribution from the 
reflecting layer itself (/ x [/i — /q]), where / can be in¬ 
terpreted as the fraction of area covered by the reflecting 
layer, for a pure broken cloud model, or as the reflectiv¬ 
ity of the cloud, for models approximating translucent 
clouds. Successive lines have a similar structure, except 
that the atmospheric contribution is only from scatter¬ 
ing between a layer and the next highest layer, and the 
contributions from both cloud and atmosphere are at¬ 
tenuated by a transmission factor TiT 2 ..Ti. In the pure 
broken cloud model this is just the fraction that is not 
blocked by the opaque elements in the previous layers, 
i.e. (1 — /i)(l — / 2 )--(l — fi)- In the translucent ap¬ 
proximation we use Ti Ri exp[—r'(l//x -|- l/^o)]; where /x 
and fiQ are cosines of observer and solar zenith angles, 
r' = (1 — + gi))Ti, Qi is the scattering asymmetry 

parameter for the ith layer, and Wi is the corresponding 
single-scattering albedo. 

The formula for the effective attenuation optical depth 
for each layer (r') is an empi rical result optimized for low 
opacity layers considered bv iSromovskv and FrvI (1200711 . 
Here, we make the simplifying assumption that the over- 
lying layers are conservative {w = 1) and symmetric scat- 
terers (<7 = 0), so that = Ti l"!. Under the same approx - 
imation we have ~ 4 x (ISromovskv and (200711 . 
so that t[ rs 2 X /i. Because our analysis of each spec¬ 
trum is at a fixed viewing geometry, we ignore possible 
angular variation in effective cloud fractions. However, 
we do need to parameterize variations with wavelength, 
which is accomplished by setting fi = /o x (A/Ao)"S 
where large particle clouds would likely have exponents 
near zero, and very small particles would likely approach 
the Rayleigh limit of = — 4. 


3.3.2. Reflecting layer I/F computations 

A comparison of single-scattering and multiple¬ 
scattering reflecting layer calculations shows that atmo¬ 
spheric multiple scattering is insignificant over the range 
1.1 < A < 4 /im, and even the 1.05 gm window adds 
no more than a few percent to the I/F obtained from 
single-scattering. The computed clear-atmosphere I/F 
values throughout this range are well below the observed 
values by 1-3 orders of magnitude, and thus particulate 
scattering is the dominant influence at all wavelengths 
considered in this analysis. For our reflecting layer cal¬ 
culations, atmospheric scattering plays a minor role that 
is adequately modeled with single scattering. 

We computed reflecting layer I/F values for unit- 
albedo and zero-albedo reflecting layers placed at 91 pres¬ 
sures logarithmically distributed between 0.5 mb and 40 
bars, using the single-scattering equations 


1 'K 

7“ = ^ ^ WiPR{es) exp(-mTi_i) 

^ i=l 

x(l — exp[—mATi])/(m/x) 

+aexY){-mTi)pL^~^ ( 2 ) 

ATi=Ti -Ti-l (3) 

— (^Ti,scat 'fi—l^scat) / (4) 

m = (l//x-b I//X 0 ) (5) 


Pr{(^s) 


3 [1 -I- J -I- (1 - (5) cos^(0s)] 
2 2 + S 


( 6 ) 


cos( 6 's) = cos((/)) v^l - - /xg - /x/io (7) 

where /x and 7 x 0 are observer and solar zenith angle 
cosines, (j) is the azimuth angle between incident and 
scattered directions, Wi is the effective single-scattering 
albedo of the gas between pressures Pi and Pi-i, a de¬ 
notes normal albedo of the reflecting layer (both 0 and 
1 are used), K defines the limb darkening character 
of the layer (assumed to have the Minnaert form in 
which reflected radiance at zenith cosine u is propor¬ 
tional to incident irradiance times Pr{9s) is 

the scalar Rayleigh phase function for anisotropic ran¬ 
domly oriented molecules, evaluate d at depolarization 
factor S = 0.020 (|Sromovskvll2005bll and scattering an¬ 
gle 0 s, Ti is the vertical extinction optical depth from 
the top of the atmosphere to pressure pi, and Ti^scatt is 
the corresponding scattering optical depth. Our calcula¬ 
tions cover the spectral range from 1.05 /xm to 3.3 /xm 
with a uniform step size of 5 cm“^, the sampling inter¬ 
val of the CH4 absorption model. For wavelengths where 
both methane and ammonia absorptions are important, 
opacities are computed 100 times for each wavelength 
(for each of the 10 terms of the correlated-k model of 
CH4, the calculation is done for each of the 10 terms 
of the correlated-k model of NH3). The most accurate 
I/F is obtained by solving the radiation transfer prob¬ 
lem for each opacity combination and then computing 
the weighted sum of the 10 or 100 calculations per wave¬ 
length. 

3.4. Multiple scattering methods 

For multiple scattering calculations, which are needed 
to accurately mod e l tran slucent clouds, we followed 
ISromovskv and Fr^ (120101) . In br ief, we used the d ou- 
bling and adding code described bv iSromovskvl (l2005ai rhh. 
We used a grid of 44 pressure levels from 0.5 mb to 10 
bars, distributed roughly in equal log increments, except 
with finer spacing in the ammonia condensation region. 
We generally used 10 zenith angle quadrature points per 
hemisphere (NQUAD) and handled the sharp forward 
scatter ing peak of large r particles using the d-Fit proce¬ 
dure of lHu et ^ (j200f)l l. The backscatter phase function 
of larger particles is not well characterized with NQUAD 
of 10 , so that calculations for low phase angles are sub¬ 
ject to errors, as shown in Section 15.11 However, it is 
unlikely that these larger particles are spherical in any 
case, so that a truncation of the spherical backscatter 
peak may actually be more realistic than keeping it. Re¬ 
duced backscatter (re lative to sph eres) is characteristic 
of aggregate particles iWestI (llQQlh. and many r andomly 
ordered crystals and spheroids (|Yang et al.l2(30Cill. includ¬ 
ing th e prolate spheroidal shapes favored bv lWong et al .1 

(Hooi. 

The most significa nt difference between the a pproach 
we used and that of ISromovskv and FrvI (j20inl l is that 
in modeling VIMS spectra we made an additional ap¬ 
proximation to account for the VIMS spectral resolu¬ 
tion being less than our model resolution. We approx¬ 
imate the line-spread function of the VIMS instrument 
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as a Gaussian of FWHM that is wavelength dependent 
(Kevin Baines, private communication 2010), and is typ¬ 
ically about 0.015 p,m. We then collect all the opacity 
values within ±FWHM of a VIMS sample wavelength, 
weight those according to the relative amplitude of the 
line-spread function, then sort and refit to ten terms in 
the exponential sums representing transmission. 

3.5. Fitting cloud models to observations 

Our cloud models were constructed as an assem¬ 
blage of discrete compact layers. Model parame¬ 
ters were adjusted to minimize using a form of 
the Levenb e rg-Ma rquardt algorithm, as described by 
I Press et al.l (Il992f) . The parameter uncertainties ob¬ 
tained from this method are based on the assumption 
of normally distributed errors, and thus cannot be abso¬ 
lutely relied on. For the reflecting layer models, we have 
supplemented these estimates by varying each fitted pa¬ 
rameter about its best fit value, refitting all the other 
parameters to minimize then setting the parameter 
uncertainty equal to the deviation needed to increase 
by one unit. This alternate method generally agreed with 
that obtained from the Levenberg-Marquardt algorithm 
to within a factor of two. It was not practical to use this 
alternate method for multiple scattering fits because of 
the enormously increased computation time that would 
be required. 

In computing y^ we tried to account for measurement 
errors as well as errors in the modeling of atmospheric 
opacity. For the measurements we assigned a random er¬ 
ror for all sources a value of 2.5% of the VIMS I/F value, 
with an additional I/F offset error of 5 x 10“^, as dis¬ 
cussed in Sec. 12.51 A more important error source, which 
is not strictly random, but can vary from wavelength 
to wavelength, is the uncertainty in radiation transfer 
modeling due to the uncertainty in opacity calculations. 
The transmission curves on which the methane opacity 
mode ls are based have an assumed fractional error of 
0.05 (|Irwin et al.ll200(il l. To convert this uncertainty to 
an I/F uncertainty, we assume that the derived optical 
depths themselves have a fractional error of similar size. 
We assigned a value of 6 % to the error in computed opti¬ 
cal depth, but converted this to I/F under the crude as¬ 
sumption that I/F would be proportional to the negative 
exponential of optical depth. This leads to a fractional 
I/F error estimate of 1-(I/F)“, where a is the fractional 
error in optical depth. Thus a 6 % error in optical depth 
would lead to errors that range from nearly zero at low 
absorption optical depths to 6 % near r = I and to 100% 
for T I. 

To judge whether one fit is significantly better than 
another on the basis of a difference in y^, we need to 
know how that difference compares to the uncertainty in 
y^ itself. The latter uncertainty arises from measurement 
and gas opacity modeling errors, and is easily determined 
from a simple Monte-Carlo calculation. For our reflect¬ 
ing layer models, we fit 94 measurements using 10 fitted 
parameters (typical), which implies an expected y^ of 
84 and an expected standard deviation of (T^2=13. The 
expected values for the multiple scattering fits are in¬ 
creased to y^ = 110 and cr^2=15 (more spectral samples 
are fit). Deviations from the expected value of y^ are a 
result of a combination of factors: statistical variability 


in y^ itself, errors in estimating uncertainties, and errors 
in the physical features of the model. 

4. EVIDENCE FOR JUPITER’S 3-/iM ABSORPTION 

To clearly display the need for a 3-/im absorber, we first 
fit VIMS spectra with aerosol models that did not con¬ 
tain such an absorber. We chose the low-latitude VIMS 
spectra from pixel locations given in Fig. [2] (low phase 
angles) and in Fig. [3] (intermediate phase angles). The 
model spectra are compared with observation in Figs. [S] 
and [TOl respectively. In these figures the observed spec¬ 
tra are shown in black, with models shown in gray. The 
three spectral regions used to constrain the fits are in¬ 
dicated by gray bars plotted in the bottom panels. The 
region from 2.85 to 3.3 fain is omitted from the fits be¬ 
cause it contains local absorption and including it greatly 
worsens the fit in other regions. The region near 2 /iin 
is avoided because of a potential connection between ab¬ 
sorption at 2 fiia and absorption near 3 /im if the cloud 
absorber is ammonia ice. The 1.60-1.68 /rm region is also 
excluded because of the previously noted VIMS respon- 
sivity correction error in this region. The fits outside 
these excluded regions are generally excellent, with y^ 
values close to expected values. 

Each layer in our 3-layer model has three parameters: 
a pressure p, a fractional reflectivity f at 2 pm and a 
wavelength dependence exponent, n. The fit parameters 
are listed in Table [T] for each of the spectra, and the 
listed uncertainties show that the parameters are well 
constrained by the observations. The need for significant 
absorption in the 2.85-3.2 pm region is apparent from the 
large ratios between model and observed spectra in this 
region (the ratios between the gray and black curves in 
Figs, ini and ITUl) . The quantitative ratios at 2.04 pm, 
2.96 pm, and 3.0-3.1 pm are listed in Tabled We see 
that as clouds brighten (I/F at 1.96 pm increases), the 
absorption needs to increase in the 3-pm region, while 
the fractional absorption near 2 ^m decreases slightly, 
suggesting that the two absorptions are not both due to 
the same cloud component (or same material). What 
appears to be happening here is that as the absorbing 
cloud layer thickens, it gets brighter in spectral regions 
where it doesn’t absorb much, but it can’t brighten in 
regions where it absorbs strongly (there its reflectivity 
saturates at a very low optical depth). Also, we find 
that the peak in the 2-pm absorption is nearer to 2.04 
pm than it is to the NH 3 peak at 2.0 pm, and the local dip 
in reflectivity near 2.97 pm is shifted from the minimum 
single-scattering albedo of pure NH 3 , which is located 
at 2.95 /xm(the peak in the imaginary index is actually 
located at 2.965 pm). 

The overall character of the 3-layer model fit param¬ 
eters is illustrated in Fig. [TTJ A high altitude aerosol 
is needed to account for the I/F in the 2.3-pm region, 
but there are two problems with the result. The effec¬ 
tive pressure of this layer is only well constrained for the 
IIxI4y, I2xI6y, and 28x3Iy spectra, for which P ^50 
mb. At that level the reflectivity is ^0.008 (r ^0.03 
for conservative symmetric scatterers), which is nearly 
20 time s the integrated value obtained by iBanfield et al.l 
(|I998al) . Banfield et al. measured low-latitude I/F val¬ 
ues of 2-3x10“^ at 2.3 ^m, while the VIMS values in 
Fig. [3 are ^10 times larger and are likely due in part 
to a VIMS measurement artifact, possibly due to scat- 
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Fig. 9.— Selected low phase angle VIMS spectra (black) compared to reflecting layer model spectra (gray). Pairs of lines surrounding 
the model spectra indicate uncertainties in opacity calculations; measurement uncertainties (not shown for clarity) are similar but slightly 
smaller. The gray horizontal bars in the bottom panel indicate the ranges over which observations were used to constrain the models. Each 
panel legend provides pixel coordinates in the VIMS cubes (as in Fig. n and values (80±13 is expected). Vertical dotted lines indicate 
pure NH 3 ice absorption features at 2.0 and 2.95 fim. 


tering of light from the VIMS grating, as described in 
Sec. 12.41 Our previous comparison with NICMOS obser¬ 
vations (Fig. [7]) also indicates that a simple offset error 
is not consistent with the observations. The substantial 
difference between the iBanfield et al.l (jl998aD reflectivity 
and our value in the 400-mb layer might also be a result 
of excessive VIMS I/F values at low levels of I/F. 


The reflecting layer fits put the second layer in the 340- 
600 mb region, with a reflectivity of 15-21%. The third 
layer, located in the 770-1600 mb region has a relatively 
high reflectivity of 40-60%. The best-fit wavelength de¬ 
pendence exponents are roughly in accord with expecta¬ 
tions. The high altitude contributions have an exponent 
of -2 to -2.9, indicating relatively small particles. This 
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As in Fig. [9] except that the spectra are for medium phase angles (as in Fig. Ell. 


Fig. 10. 

might have been closer to -4 if the excess I/F near 2.3 
/im could have been removed. The next layer has best-fit 
exponents in the -0.75 to -1 range, which indicates larger 
particles, and the deepest layer has even smaller expo¬ 
nents (-0.5 to -1-0.15), indicating even larger particles. 


We also tried fitting the observations with gas opacity 
profiles other than the nominal EZ profile, but even for 
the least opaque region sampled in the NEB, the best fit 
was obtained with the EZ profile. 

The latitude dependence of the fit parameters (Eig. 
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TABLE 1 

Best fit parameter values eor 3-layer reelecting-layer 
EITS TO SELECTED VIMS SPECTRA. 


Spectrum 

Parameter 

Layer 1 

Layer 2 

Layer 3 

llxl4y 

p (bars) 

"VoisoT^^ 

0 3Q9+o-u:i! 

U.Oy^_Q 

0-919l^;“5 

llxlSy 

p (bars) 

0 .0020l°;°o? 

0.380l°;“ 

0.93Sl°;°° 

12 xl6y 

p (bars) 

o.osastnis 

0-4iVS:o3 

0.976j;o)f 

13xl4y 

p (bars) 

o.ooos 

0.602tO;°) 

1.624l°;)° 

25x35y 

p (bars) 

o.ooos 

OA62t°o°l 

1.023l°;)° 

26x31y 

p (bars) 

0.0179tni? 

0.397t°;“ 

0.7S4l°;°° 

27x34y 

p (bars) 

0.0019l°;°o? 

0.343l°;°l 

0.812l°;J° 

28x31y 

p (bars) 

0.05S8tni4 

n q'77+0.03 
'-0.03 

n 779+0.09 

y- > <^_o.08 

llxl4y 

/ 

0.0070t°;°°22 

0.183l°;“( 

0.6241°;°° 

llxlSy 

/ 

0.00S8t°;°°°® 

0.1771°;“ 

O.S9ll°;“ 

12 xl6y 

/ 

0.0082t°;°°® 

0.146l°;°) 

0.3681°;“ 

13xl4y 

/ 

U.UUOU_Q QQQ2 

0.209l°;“ 

o.eoilg'i^ 

25x35y 

/ 

U.UUZ^_Q QQQ2 

0.138l°;°( 

0.3731°;“ 

26x31y 

/ 

U.UUZO_Q QQQQ 

0.184l°;“ 

0.2481°;“ 

27x34y 

/ 

0 0037+0 0019 

U.UUOf _0.0011 

0 -IQO+O OO 

u.iyz_Q 00 

0.4221°;“ 

28x31y 

/ 

o 

b 

o 

CD 

o 

1 + 
o o 

b b 

o o 

o o 

U.Zi4_o 01 

0.46Sl°;“ 

llxl4y 

n (in A") 

-2.86t°oit 

0 QO + O-IO 

G-y-^-O.ig 

0 1 3+0-23 
-U.10_o.i8 

llxlSy 

n (in A") 

-2.06l“;33« 

i.UU_o 21 

^•1^-0.10 

12 xl6y 

n (in A") 

-L-»^_0.44 

-1-1+ofo 

-0 24’*''^'^^ 
^•^‘*_0.39 

13xl4y 

n (in A") 

2 34+0-12 

-0.6Sl°;“ 

n 31+0.57 
-U.Oi_o.55 

25x3Sy 

n (in A") 


-o.9il°:)« 

0 26+*^'^^ 
-U.ZO_o.38 

26x31y 

n (in A") 

-2.86l°;g 

- 0 . 8 sl°;)? 

0 11+0.08 
-U.li_o.o8 

27x34y 

n (in A") 

Q -1 y-|-0.29 
_0.30 

0 oc + O.lO 
-u.oo_o.io 

0 1 ir + 0.39 
^•10-0.47 

28x31y 

n (in A") 

-2 34+0-25 

0 yir + o.o? 
y- <3_0.07 

0 03+0.59 
G.UO_o 55 


NOTE: Layer-1 pressures without uncertainties had best fit val¬ 
ues at the lower limit of the fit range. 


[m displays an increase in pressure towards the north 
equatorial belt for both layers 2 and 3. (The Cassini 
VIMS 0.45-/tni image in Fig.[5p indicates that the NEB 
extended from about 8 °N to 18°N in 2000.) The trends 
in Layer 1 pressures is not shown because it is likely not 
meaningful (several fits yielded best-fit pressures at the 
lower limit of the allowed range). The pressure changes 
seem to have more impact on the observed I/F than the 
changes in optical depths. 

The pressure level at which the 3-^m absorber resides 
is roughly that of layer 2 in the 3-layer model, which is 
in the 300-600 mb range where NH 3 ice is a plausible 
condensate. The location is inferred from the cloud re¬ 
flectivity and the penetration depth profile given in Fig. 
[TJ Vertical penetration depths at 2 /rm and 3 /rm are 
comparable and limited to pressures less than 1 bar or so. 
The relatively high reflectivity seen at 2 /xm (compared to 
that at 3 /xm) is produced by contributions mainly from 



Fig. 11. — Pressure and reflectivity parameters for 3-layer reflect¬ 
ing layer models that provide best fits to VIMS spectra, as given 
in Table ^ The dash ed curve displays the vertically integrated 
reflectivity derived by IBanfield et al.l ijlQQSal i for the GRS but is 
broadly representative of other regions as well. 




Fig. 12.— Pressure and reflectivity parameters versus latitude for 
best-fit 3-layer reflecting layer models given in Tabled Parameters 
are plotted for layer 1 (diamonds), layer 2 (open squares), and layer 
3 (filled circles). 

layer 2, which would be seen even better at 3 /rm. It 
thus follows that layer 2 must be strongly absorbing at 3 
fim to account for its much lower I/F at that wavelength. 
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TABLE 2 

Fitted spectra compared to VIMS spectra in regions of enhanced absorption. 


Spectrum 

Measured 

I/F at 1.95 fim 

Ratio of model I/F to measured I/F 
at 2.04 /rm at 2.96 /rm at 3.0-3.1 ^m 

13xl4y 

0.048 

1.300 

3.44 

2.76 

25x35y 

0.058 

1.261 

4.32 

3.35 

12 xl6y 

0.095 

1.191 

4.03 

3.40 

llxl5y 

0.124 

1.233 

5.45 

4.24 

llxl4y 

0.125 

1.211 

6.65 

4.68 

26x31y 

0.133 

1.200 

5.43 

4.06 

27x34y 

0.137 

1.180 

4.55 

4.36 

28x31y 

0.163 

1.143 

6.54 

5.02 


Note: above rows are in order of increasing I/F at 1.95 fim. 


However, the layer-2 and layer-3 fractions in this model 
are too large for the reflecting layer equations to provide 
physically meaningful parameter values if the clouds are 
actually uniform and translucent rather than broken and 
opaque. For the case of uniform translucent clouds, we 
need multiple scattering fits, which are presented after 
first identifying a plausible composition for the particles 
in such clouds. 

5. MODELING THE S-pm ABSORPTION. 

5.1. Candidate 3-^m absorbers. 

Candidate cloud materials expected in Jupiter’s at¬ 
mosphere that also absorb light in the 3-/im region in¬ 
clude NH 3 , NH 4 SH (ammonium hydrosulfide), N 2 H 4 
(hydrazine), and water ice. Their real and imaginary 
refractive indexes are shown vs. wavelength in Fig. 1131 
Because of its very high imaginary index, water ice can¬ 
not produce the absolute reflectivity level needed at 3 pm 
and the transparency needed at 2.7 pm. Trial fits with 
water ice resulted in very large values. Water is also 
a poor candidate because of its extremely low mixing ra¬ 
tio in the pressure range where the 3-pm absorber seems 
to be located. Hydrazine has a pair of strong absorption 
peaks between 3 and 3.2 pm, which were very apparent in 
in preliminary model spectra, and thus seemed inconsis¬ 
tent with observations. Among the more plausible com¬ 
pounds, only NH 3 has a known absorption at 2 pm (the 
2.25 pm absorption feature of NH 3 would not be visible 
due to overlying gas absorption). NH 4 SH has an absorp¬ 
tion that is roughly comparable to that of NH 3 at 3 pm 
(except for the sharp NH 3 feature at 2.95 pm), but con¬ 
tinues to increase beyond 3.1 pm where NH 3 absorption 
drops significantly. In the following detailed discussion 
we only include ammoni a and ammonium hydrosu lfide 
as viable candidates (see iSromovskv and Fr^ (1201011 for 
an evaluation of a broader range of materials). 

The cloud layers at which the 3-pm and 2 -pm absorp¬ 
tions contribute do so primarily through their reflected 
intensities. Thus, it is useful to consider the spectral re¬ 
flectivity of unit optical depth layers of NH 3 and NH 4 SH 
for a variety of particle sizes and for low and medium 
phase angle observations. As shown in Fig.[T31 the reflec¬ 
tivity plots provide clues as to what sorts of layers can 
provide the appropriate shortwave reflectivity and the 


needed 3-pm absorption. Fig. [T3] also displays the effect 
of NQUAD (quadrature angles per hemisphere) on the 
accuracy the computed reflectivity. The difference be¬ 
tween calculations for NQUAD values of 10 and 60 show 
that as particle size increases much beyond a few microns 
model computations for low phase angles are much more 
demanding than for medium phase angles. This occurs 
because the scattering phase function is very smooth at 
middle phase angles (and scattering angles) but becomes 
quite sharp and more strongly wavelength dependent for 
large size parameters in near backscatter configurations 
(low phase angles). 

A prominent feature in the reflectivity of ammonia lay¬ 
ers is the deep and sharp minimum at 2.95 pm. The lack 
of a such a feature in the VIMS spectra, or in the hi gher 
resolution ISO spectrum (ISromovskv and 1^1201011 im¬ 
plies that the top cloud layer on Jupiter cannot be made 
of pure NH 3 ice. If only ammonia is considered for the 
middle layer (500-700 mb), then only relatively large par¬ 
ticles provide roughly the correct spectral variation in 
reflectivity. But even then, to fill in the deep minimum 
at 2.95 pm some additional reflectivity must be provided 
by a gray particle layer extending above the ammonia 
layer. Small particle layers of ammonia (0.25-1 pm in ra¬ 
dius) are relatively too bright just beyond 3 pm, and they 
also contribute too much reflection at short wavelengths. 
Large ammonia particles produce a relatively strong 2- 
pm feature that is too strong to be compa t ible w ith ob¬ 
servations, as pointed out by lirwin et all (|2001h . The 
NH 4 SH layers provide a flatter absorption between 3 and 
3.4 pm, with a slope between 3 and 3.1 pm that depends 
on particle size. Note that NH 4 SH particle layers are 
brighter than NH 3 layers of the same optical depth and 
particle size, a result of the higher real index of NH 4 SH. 

The VIMS spectra display a significant variation in 
brightness at 1.95 pm, without a significant variation in 
brightness at 3 pm, and the fractional absorption at 2 
pm remains relatively constant, while the fractional ab¬ 
sorption at 3 pm changes by a factor of two. Another 
feature of the spectra is the small dip at 2.97 pm, which 
is close to the 2.95-pm dip that is characteristic of an 
ammonia ice layer. These facts suggest that both ab¬ 
sorbers are present at pressure levels that are visible at 2 
pm and 2.9-3.3 pm, with NH 4 SH providing most of the 
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Fig. 13. — Re al (bottom) and imaginary (top) components of the refrac tive index vs. wavel ength for cand i date 3 -m absorbers, including 
H 2 O results of IWarr^ ll 198411 . NH 3 results of IMartonchi^_gt_^lJJ[198^ , NH 4 SH results of IHowett et al.1 II2007I 1 and Barbara Carlson 
(personal communication), and N 2 H 4 results of lcTaDp'^inT^illerl 11199111 1. 


cloud mass and most of the 3-/im absorption, and ammo¬ 
nia ice contributing a small amount of 2 /rm absorption 
and helping to produce the small dip at 2.97 /im. These 
suggestions are quantitatively evaluated in the following 
section of multiple scattering model fits. 

5.2. Models with an ammonia cloud layer. 

We first considered the possibility that NH 3 could pro¬ 
vide all the needed 3-/im absorption. Our initial multi¬ 
ple scattering models consisted of four scattering layers. 
The top layer was a small particle haze characterized by 
symmetric Henyey-Greenstein (H-G) phase function, a 
wavelength dependence, an adjustable optical depth, 
and a fixed pressure of either 15 mb or 50 mb, with the 
choice guided by the reflecting layer fits. This top layer is 
needed to provide reflectivity in the 2.3-/rm region. The 
second layer is a Mie-scattering layer of conservative par¬ 
ticles with adjustable particle size, pressure and optical 
depth. The third layer consists of spherical NH 3 ice par¬ 
ticles, with adjustable particle size, pressure and optical 
depth. The fourth layer was modeled as a nearly conser- 
vative layer {w = 0. 9 97) u sing; the H-G phase function of 
iSromovskv and Fr^ (I2002D and a A" wavelength depen¬ 
dence, with fit-adjustable pressure, optical depth, and A 
exponent n. We also found comparable fits treating this 
as a Mie-scattering layer of NH 4 SH particles, although 
in that case a stronger backscatter peak is needed to 
provide enough reflectivity. Without this deeper layer 


it appears impossible to match the spectral shape of the 
reflectivity peak near 1.23 /rm. This structure, with NH 3 
as the only aerosol absorber, provides an excellent fit to 
the lower opacity regions, such as the 13xl4y spectrum, 
as shown in Fig. I151 A.B. In this example the ammonia 
layer pressure (673 mb) is high enough that its spectral 
features at 9.6 yrm and 26 /rm would likely be hidden 
from view in observed thermal emission spectra, where 
such features are not easily observed. 

In regions of higher opacity, such as location IIxI4y 
in the low phase angle cube, the spectra cannot be as 
accurately fit with this structure. As shown in Fig. [T5G 
and D, there is an absorption at 2 /im in the model spec¬ 
tra that does not appear in the observed spectra. The 
best-fit parameter values for both fits are listed in Table 

El 

5.3. Models with NH 4 SH as the only 3-^m absorber. 

We next considered the possibility that NH 4 SH could 
provide all the absorption at 3 /im. This would certainly 
eliminate the problem of excessive absorption at 2 /im, 
although it would not be able to fit the feature at 2.97 
/im. For this study we used a model similar to that 
used for the NH 3 analysis, with the third layer consist¬ 
ing of NH 4 SH instead of NH 3 . We used a composite re- 
fractive index functi on for NH 4 SH: at wavelengths where 
iHowett et al.l (|2007ll observations exist we used the av¬ 
erage of their 80 K and 160 K measurements; elsewhere 
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Fig. 14.— I/F vs. wavelength for unit optical depth layers of spherical NH 3 (left) and NH 4 SH (right) particles of radii from 0.3 fim to 
10 fira. Results for two phase angles are shown: 2.5° (top panels) and 67.8° (bottom panels) and for NQUAD values of 10 (solid) and 60 
(dot-dash). Note that the NH 4 S H calculations betw een 1 and 2.5 /rm are based on interpolation of refractive index measurements at 1 and 
2.5 fim. Fhost reflection spectra IIFanale et al.|[T977ll show broad absorption features that may be significant between 1 and 2 fim, but no 
refractive index values are available. 


we used Barbara Carlson’s compilation (personal com¬ 
munication 1994). Again we tried the two extreme cases 
(13xl4y for low cloud opacity, and llxl4y for high cloud 
opacity). The best-fit spectra are displayed in Fig. [1^] 
and the corresponding parameters given in Table [3] The 
low-opacity fit is slightly worse than the NHa-only fit 
(by only 0.26 x 0 -^ 2 ), but the high opacity fit is better (by 
lx 0 -^ 2 ) than we obtained using NH 3 as the sole 3-/im 
absorber. However, the NH 4 SH fits have specific local 
deficiencies: as expected, they don’t fit the feature at 
2.97 fim, they still don’t provide close fits in the 1.9-2.0 
fim region, and they have problems fitting the 2.7-fim 
peak and its long-wavelength shoulder. 

There is also some question concerning the appropriate 
NH 4 SH absorption model to use. The standard equilib¬ 
rium chemistry model places NH 4 SH condensation near 
160 K, and particles formed at that level and lofted to 
higher altitudes will retain the same crysta l structure, as 
eviden t from laboratory measurements of iFerraro et al.l 
(|1980r) . On the other hand if particles formed at some¬ 
what lower temperatures, a different crystal structure 
would apply and the absorption would be close to that 
measured at 80 K. By picking an average of measure¬ 
ments at 80 K and 160 K, we are allowing for a 50:50 
mix of formation temperatures above and below 160K. 
A trial spectral fit to the llxl4y spectrum, assuming that 
all the NH 4 SH was formed at high temperatures and that 
NH 4 SH was the sole absorber, was somewhat degraded: 
increased by l.5xa^2. A trial spectral fit using the 


80 K absorption model for NH 4 SH, was also worse, but 
less significantly so (x^ increased by 0 . 8 xo-^ 2 ). Thus the 
observations slightly favor a mix of formation tempera¬ 
tures. 

We see that neither NH 3 alone nor NH 4 SH alone pro¬ 
vides accurate fits to spectra over the full range of cloud 
structures. NH 4 SH seems useful for fitting the high opac¬ 
ity cases and either NH 3 or NH 4 SH are equally use¬ 
ful for fitting the low-opacity cases (the improvement 
of NH 3 over NH 4 SH in those cases is really insignifi¬ 
cant). It is worth noting however, that NH 4 SH and NH 3 
have somewhat complementary fitting defects: (1) NH 3 
tends to be low near 2.85 fim where NH 4 SH tends to 
be slightly higher (see location Y in Figs. [15] and 1161) : 
(2) NH 3 tends to be increasingly high beyond 3.05 fim 
where NH 4 SH tends to be increasingly low. A combi¬ 
nation of the two absorbers thus offers the possibility of 
reducing the net fitting defects in these areas. But the 
most compelling reason to keep both NH 4 SH and NH 3 
particles in the model is that when NH 3 is combined with 
a non-absorbing material it becomes possible to fit the 
2.97-/rm feature (location Z in Figs. and ITBl). While 
the NH 3 reflectivity minimum for pure NH 3 particles is 
displaced from this feature, a shift towards the peak in 
the NH 3 imaginary index (at 2.965 fim) can be accom¬ 
plished if NH 3 is diluted by coating or mixing it with a 
conservative material, or by condensing it as a coating 
on a conservative core, as will be shown in the following 
section. 
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1.0 1.5 2.0 2.5 3.0 


Wavelength, /xm 


Fig. 15.— Multiple-scattering 4-layer fits to lowest (A) and highest opacity (B) cloud structures at low phase angles. These fits use 
NH 3 Mie particles for the third layer, which fits well only for low opacity clouds. Model spectral points are displayed as filled circles, 
measurements as solid lines with gray lines indicating uncertainties. Vertical dotted lines at 2.0 and 2.95 fim mark locations of minima 
in the single scattering albedo NH 3 ice for particle radii from 0.2 to 20 fim. Companion plots show the difference between model and 
measurement divided by the expected combined error. Letters X, Y, and Z mark locations of noteworthy differences. 


5.4. Models with both NH 3 and NHiSH layers. 

We first considered a structure similar to the previous 
models except that we inserted a layer of NH 3 parti¬ 
cles between a conservative Mie layer and a Mie layer 
of NH 4 SH pa rticles, which is esse n tially the same struc¬ 
ture used by iSromovskv and Fi^ (|2010ll to successfully 
fit the ISO spectrum. This approach produced relatively 
good fits for the low-opacity spectra, such as 13xl4y and 
25x35y, with best fits obtained with all three particle 
layers at nearly the same pressure. However, this struc¬ 
ture did poorly at fitting the high-opacity spectra, such 
as llxl4y, for which the result was a poor match to the 
2.97-/im feature and too low an I/F in the 1.95-2.1 /rm 
region. This suggests that the structure inferred from 
the ISO spectrum does indeed not apply at all locations, 
which is not surprising, given the large FOV of the ISO 
spectrum and the wide variety of clouds averaged to¬ 


gether. 

We then consi d ered the possibility raised by 
ISromovskv and Fi^ (|2010ll that an NHa-coated particle 
might provide a better fit, either as a coating on the 
larger NH 4 SH particles, or as a coating on the smaller 
conservative particles usually found near 350 mb. We 
computed scattering properties of coated spheres using 
the Fortran subroutine DMiLa y, based on original work 
by iToon and Ackerniaiil (1198111 , and further refined by 
W. Wiscombe (program code and documentation avail¬ 
able at ftp://climatel.gsfc.nasa.gov/wiscombe/). Al¬ 
though the 2.95-/rm absorption feature of pure NH 3 is 
not seen in the VIMS spectra, the feature near 2.97 /rm 
could still be due to NH 3 when it is applied as a coat¬ 
ing. An example is provided in Fig. 1171 The upper right 
panel of that figure displays the reflectivity of a thin layer 
of spherical particles of 0.8 fj.m in radius, viewed in a 
backscatter configuration. When the particles are com- 
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TABLE 3 

Best-fit cloud structures with single and multiple 3-/jm absorbers. 


Parameter 

Low Opacity Spectrum (13xl4y) 

High Opacity Spectrum (llxl4y) 

Layer-2 Composition 

1.40 -1- Oi 

1.40 -1- Oi 

NH3 shell* 

1.40 -1- Oi 

1.40 -1- Oi 

NH3 shell* 

Layer-3 Composition 

NH3 

NH4SH 

NH4SH 

NH3 

NH4SH 

NH4SH 

Pi (bars) 

0.015 

0.015 

0.015 

0.050 

0.050 

0.050 

P2 (bars) 

0.517 

0.553 

0.546 

0.303 

0.307 

0.336 

P3 (bars) 

0.628 

0.608 

0.566 

0.625 

0.646 

0.635 

P4 (bars) 

1.261 

1.158 

1.128 

1.616 

1.285 

1.252 

n 

0.014 

0.014 

0.014 

0.023 

0.023 

0.025 

T2 

0.270 

0.369 

0.414 

0.306 

0.291 

0.420 


1.569 

0.386 

0.231 

6.63 

2.83 

2.10 

n 

4.66 

5.51 

5.19 

14.50 

20.0 

10.64 

Layer-1 A exponent 

-3 

-3 

-3 

-3 

-3 

-3 

r2 (/rm) 

0.289 

0.293 

0.303 

0.265 

0.237 

0.259 

rs (/rm) 

6.98 

13.35 

6.86 

6.51 

12.54 

10.80 

Layer-4 A exponent 

-1.75 

-1.21 

-1.25 

-1.28 

-0.625 

-1.12 


125.79 

129.70 

118.06 

175.41 

160.47 

124.44 


Notes: For layer 1 we used a symmetric phase function and fi xed pressures and A exponen ts. For layer 
4 we used the double Henyey-Greenstein phase function of ISromovskv and FrvI 11200211 . Parameter 
uncertainties are given for the composite fit in Table |4] Optical depths are evaluated at 2 pm. The 
NH3 shell* is a simulation of an NHa-coated particle as described in the text, differences less 
than 15 are of marginal significance. 


posed of pure NH 3 , a deep minimum is seen in the layer 
I/F, exactly at 2.95 /xm (shown as the solid curve). But 
when the particles are composed of N=1.75 material out 
to 80% of the radius, and the remainder consists of NH 3 , 
the I/F minimum is shallow and shifted close to the 2.97- 
/xm peak in the imaginary index of NH 3 . A similar effect 
is seen when an ammonia core is covered by a heavy coat 
of benze ne (thickness equal to co re radius), as shown in 
Fig. 5 of lKaloeerakis et all (I2008D . It is also noteworthy 
that the coated particle has a greatly reduced extinc¬ 
tion efficiency relative to pure NH 3 particles at thermal 
wavelengths, making them less likely to be detected in 
thermal spectra. 

Also shown in Fig. [T7]is the I/F for uniform particles 
with a synthetic refractive index that roughly simulates 
the reflectivity of the composite particle. This allows 
us to incorporate the essence of a coated particle into 
the Mie calculations contained in our multiple scatter¬ 
ing model without the time penalty of carrying out the 
more complex and time-consuming coated sphere calcu¬ 
lations. We created the synthetic imaginary index by in¬ 
serting a Gaussian absorption feature of amplitude 0.08 
and FWHM of 0.06 p,m at 2.98 fxm, and using 1/2 of the 
NH 3 imaginary index, wherever that exceeded the wings 
of the absorption feature. We then using the Kramers- 
Kronig relation, following iHowett et al.l (I2007D . to com¬ 
pute a spectrally varying real index that was consistent 
with the imaginary index. The synthetic index provides 
an absorption feature that crudely matches the coated- 
sphere feature. However, the detailed character of this 
absorption depends on particle size, core material, shell 
thickness, and probably particle shape. Thus, the ex¬ 
ample shown in Fig. [T3 should be considered suggestive, 
rather than typical, of the absorption feature that can 
be created by coatings of NH 3 on other materials. We 


tried to create absorption features of this type on larger 
(10-15 ^m) particles with NH 4 SH as the core material, 
but were unsuccessful, largely because of the rapid in¬ 
crease in absorption by NH 4 SH on the long-wavelength 
side of the feature. We thus find that the most plausible 
means of creating an absorption feature at 2.97 /xm is 
by coating small conservative particles with NH 3 (or by 
coating NH 3 cores with a conservative shell). We thus 
proceeded to carry out fits using small particles made 
of the synthetic index material that roughly simulated 
these composite particles. 

Our revised structure placed the synthetic index parti¬ 
cles in layer 2 , to simulate the absorption by an ammonia 
coating on a conservative core. Fit-adjusted parameters 
included the pressures of the bottom three layers, the 
optical depths of all four layers, and the particle sizes of 
layers 2 and 3. From the comparisons of model and mea¬ 
sured spectra, shown in Figs.fTBlandfTOl we see that these 
model fits are much better than were obtained when only 
a single 3-^m absorber was used . The shape of the 
2.7-fxm peak and its shoulder down to the 2.97-p,m NH 3 
feature is especially well matched, and overall values 
are much improved as well, as can be seen in Table [31 
which compares parameter values and fit quality for all 
three absorber models (the improvement in is 2.4-3.4 
xcr ^2 for the high opacity case and 0.5-0 .8 XtT ^2 for the 
low opacity case) . For the 13xl4y spectrum, we see 
that layer 2 and layer 3 are closer together than for the 
single absorber models, the layer 3 optical depth is less, 
and the bottom cloud pressures are very similar. For the 
llxl4y spectrum, which is from a much cloudier region, 
the pressures of layers 2 and 3 are similar to those for the 
NH 3 case, but the optical depth for layer 2 is about 40% 
more than for either of the single-absorber cases, and the 
pressure of the bottom layer has been reduced from ^ 1.6 
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As in Fig. 1151 except these 4-layer multiple scattering fits use NH 4 SH Mie particles for the third layer. 


Fig. 16.— 
bars to 1.25 bars. 

The best-fit parameter values for the two-absorber 
model are listed for all eight fitted spectra in Table 01 
We find that the preferred sizes of the NH 4 SH particles 
are relatively large (r ~ 7-16 /tm), which is comparable 
to the sizes we inferred for models in which only one ab¬ 
sorber is present. On the other hand, the preferred size 
of the NHa-coated particles is quite small, ^0.3 ^m. In 
these models the NH 4 SH layer provides most of the 3-/im 
absorption, but NH 3 provides additional spectral shap¬ 
ing that helps to achieve a better match to the observed 
spectra. Because we did not include sufficient quadra¬ 
ture points to characterize the complex backward scat¬ 
tering phase functions of the larger particles (if treated as 
spheres), the larger particle sizes are mainly constrained 
by the spectral shape of the 3-/rm absorption feature, 
rather than the large scale wavelength dependence. For 
models in which there was a close proximity of layers con¬ 
taining NH 3 and NH 4 SH particles, cross-correlations re¬ 
sulted in very poorly constrained uncertainty estimates. 
This was handled by forcing the layer 2 pressure to follow 


the layer 3 pressure by a fixed offset (usually 20-30 mb). 

The vertical structure of the two-absorber models is 
visualized in Fig. [20] by plotting layer optical depth as a 
function of pressure. Comparing this to Fig. 1111 we see 
that the derived pressures are similar, with layer 3 of the 
reflecting-layer model being split into two layers of differ¬ 
ent composition in the current model. For comparison we 
also show two independent results, converted to cumu¬ 
lative vertical optical depth. Because most of our layers 
have much higher optical depths than overlying layers, 
our cumulative optical depths would be stair-steps that 
would be just slightly to the right of our plotted layer val¬ 
ues. The pair of triple-dot-d a sh cu rves are from belt and 
zone profiles of lirwin et al.l (120011) . These are well be¬ 
low most of our values, typi cally by factors of 2 or more. 
Both our results and those of lirwin et al.l (j2001[ ) probably 
contain far too much aerosol optical depth in the strato¬ 
sphere because both NIMS and VIMS obseryations near 
2. 3 jUm seem exce s siyely h igh compared to measurements 
of iBanfield et al.l (ll998aD and to NICMOS obseryations 
(Fig. [7]). The high stratospheric opacity from our VIMS 
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2.80 2.85 2.90 2.95 3.00 3.05 
Wavelength, //m 



8.6 8.8 9.0 9.2 9.4 9.6 9.8 10.0 
Wavelength, fj,rT\ 



Fig. 17.— Comparison of scattering properties of 0.8-/rm radius spherical particles of pure NH 3 (solid curves), N=1.75 cores coated with 
NH 3 shells of 20% of total radius (dot-dash), and particles with a synthetic refractive index described in the text (dashed). Shown are 
extinction optical depth normalized to 0.34 at 2 pm (upper left), extinction efficiency near 9.4 pm (lower left) and 26 pm (lower right), 
and reflectivity of a thin layer in the backscatter direction near 2.95 pm (upper right). 


analysis is a res ult of the VIMS artifact previously dis¬ 
cussed in Sect. nn which is probably caused by light 
scattered off the VIMS grating. At least part o f the h igh 
stratospheric opacity obtained by lirwin et al.l ( 120011 ) is 
due to their assumption of a Henyey-Greenstein phase 
function with an asymmetry parameter of 0.5, which pro¬ 
vides far less backscatter than the sub-micron particles 
which are the likely source of stratospheric scattering. 
To make up for this lack of backscatter, the inferred op¬ 
tical depth for their model had to increase to match I/F 
measurements. A more detailed comparison with NIMS 
results can be found in Sec. 16.21 
The variations of model parameters with respect to lat¬ 
itude are shown in Fig. [2T] These results show some of 
the tendencies observed in the reflecting layer parame¬ 
ters. The cloud layers seem to descend to higher pres¬ 
sures moving from zone to belt, but there is not a clear 
tendency in optical depth. The remarkable feature is the 
close proximity of the layer of small NH 3 particles to the 
layer of much larger NH 4 SH particles, suggesting that 
both particle populations might actually overlap in the 
same layer, perhaps with different scale heights. The de¬ 
scent to higher pressures in the NEB also helps to hide 
the NH 3 particles from observations at thermal wave¬ 


lengths, helping to explai n the latitudinal de pendence of 
NH 3 features det ecte d bv iWong et al.l (j2004l ). as we dis¬ 
cuss later in Sec. o 

There are two counter examples to the general rule of 
near overlap of NH 3 and NH 4 SH layers. These are for two 
near equatorial spectra (llxl4y and llxl5y), for which 
these two layers are separated by nearly 300 mb. For 
these spectra the NH 3 layer stayed at low pressures char¬ 
acteristic of the equatorial zone, but the NH 4 SH layer 
moved to higher pressures, as did the bottom cloud layer. 
Optical depths for these cases did not deviate much from 
those inferred for the other spectra. 

5.5. Sensitivity of spectra to model parameters. 

The contributions of each model layer to the spectral 
features in the llxl4y spectrum are illustrated in Fig. [521 
where the best-fit spectrum is compared to four spectra 
that each have one layer removed. This shows that the 
NH 3 -coated particle layer (336 mb) provides important 
contributions in several regions, including near 1.5 /rm , 
near 2 ^m, and from 2.95 ^m to 3.2 p,m. The NH 4 SH 
layer (635 mb) provides important contributions near 2 
/rm, 2.5 /xm, and from 2.65 /xm to 2.95 /xm, but has so 
much absorption beyond 2.97 /xm that it contributes es- 
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Fig. 18. — Selected low phase angle VIMS spectra (lines) compared to translucent layer multiple scattering model spectra (filled circles) 
using both simulated NHs-coated and NH 4 SH cloud layers. Light gray lines indicate combined uncertainty bands. The VIMS spectra from 
1.58 to 1.68 fim are not used for constraining fits because of a responsivity error. Each panel legend provides pixel coordinates in the VIMS 
cubes (as in Fig. n and values (a value of 110±15 is expected). Here a vertical dotted line marks a wavelength of 2.97 / im . 

sentially nothing to the I/F from 2.97 to 3.2 /rm. Thus, mb) is seen to make significant contributions only to the 
these two different absorbers are seen to make distinctly window regions. 

different contributions. The deepest cloud layer (1252 The sensitivity of model spectra to specific model pa- 
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As in Fig. 1181 except that the spectra are for medium phase angles (as in Fig. 0. 


Fig. 19. 

rameters can be discerned from the fractional derivatives 
of a model spectrum with respect to each of the model 
parameters. This is shown in Fig. [23] for a model that is 
the best fit to the 25x35y spectrum. This example illus¬ 
trates the relative independence of the spectral deriva¬ 


tives from each other, confirming the possibility of using 
observations to independently constrain the parameters. 
For example, although the upper layer of particles is only 
important in the more opaque regions of the spectrum, 
the fractional derivatives with respect to pressure (right. 
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TABLE 4 

Best fit parameter values eor 4-layer multiple-scattering 2-absorber fits to VIMS spectra. 


Spectrum 

Parameter 

Layer 1 
(A") 

Layer 2 
(NH3 shell) 

Layer 3 
(NH4SH) 

Layer 4 
(A") 

llxl4y 

p (bars) 

0.0500 

U.OOU_Q 

0 633+0-042 
U.OOO_Q 

i-252l°;°l 

llxl5y 

p (bars) 

0.0150 

0 330+0-003 

U.OOU_Q QQ3 

u.Diy_Q 

1.268l°;)° 

12xl6y 

p (bars) 

0.0150 

0.381 (P3-0.030) 

0 411+0-006 

0.885l°;°t 

13xl4y 

p (bars) 

0.0150 

0.546 (P3-0.020) 

0.566l°°°( 

1.1281°;°) 

25x35y 

p (bars) 

0.0150 

0.460 (P3-0.060) 

0.520l°°ll 

1.007l°;°) 

26x31y 

p (bars) 

0.0500 

0.409 (P3-0.030) 

0 43Q+O-OIO 

0.954l°;°t 

27x34y 

p (bars) 

0.0150 

0.345 (P3-0.013) 

0 330+0.028 
u.ooo_Q 023 

0.8851°;°^ 

28x31y 

p (bars) 

0.0150 

0.341 (P3-0.030) 

n *^71 +0.001 
u.oii-o.ooi 

0.7891°;°? 

llxl4y 

T 

0.0250^^0! 

0.420lo°| 

^•l'J_0.49 

10.6411°^ 

llxl5y 

r 

0.0260tnol 

0.450lo°( 

^•-‘^'^-0.24 

19.991);°° 

12xl6y 

T 

0.0204t°;™2 

0.370i°;°l 

n Qq+0.02 

^•'^«J_0.02 

Q 07+0.12 

'J.ol_o.i2 

13xl4y 

r 

0.0143t°;°“ 

0-414l°°® 

Q 23+0.01 

5.l9l^;« 

25x35y 

r 

O.OlOllHo? 

0.284tO;°( 

J--^u_o.o6 

8.03lg;^^ 

26x31y 

T 

0.0168tHo2 

0-374t0;°t 

1 yr+0.24 
-*-• 'O_0.24 

6.80l°;«^ 

27x34y 

T 

0.016o1°;°09 

0.387l°;)I 

1 30+0-35 
i.c>y_o.3i 

10.90lt,« 

28x31y 

r 

0.0129tHol 

0.3781°°® 

l-48lgig 

12 1 n+O-^® 

1Z.1U_3 03 

llxl4y 

r {pm) or n 

n= -3 

r= 0.2591q Q 7 

r= 10.801°;2° 

rj— 1 12'*"*^'®^ 
n— i.iz_045 

llxl5y 

r (pm) or n 

n= -3 

r= 0.280l°;°® 

r= 16.39lli® 

n= 0.00 

12xl6y 

r (pm) or n 

n= -3 

r= 0.270l° °( 

r= 16.48l°l° 

n= -0.82l°;)? 

13xl4y 

r (pm) or n 

n= -3 

r= 0.303l°« 

r= 6.86l°;|° 

rj— 1 25"*"*^’^^ 
[L— i.zo_o 03 

25x35y 

r (pm) or n 

n= -3 

r= 0.364l°°( 

r= 11.831°;°° 

n= 0.43l°'+ 

26x31y 

r (pm) or n 

n= -3 

!■= 0.3551° °® 

r= 13.841°;°° 

n= 1.73 

27x34y 

r (pm) or n 

n= -3 

!■= 0.3331° °| 

r= 11.301°;°° 

n= 1.56 

28x31y 

r (pm) or n 

n= -3 

r= 0.304121^ 

r= 10.931°;°° 

n= 1.64l°;°° 


NOTE: Only formal uncertainties from the fits are given; see text for other sources of un¬ 
certainty. Layer 2 pressures without uncertainties follow layer 3 pressures with fixed offsets. 
Some Layer 4 wavelength exponents were also fixed to stabilize uncertainty estimates. 


top row) and optical depth (right, third row) are dis¬ 
tinctly different. A similar difference can be found for 
other layers, and especially large differences exist be¬ 
tween different layers. Note that particle size changes 
in layer 2 (the NHa-coated layer) have the greatest ef¬ 
fect at 2.97 fj,m and from there to 3.2 /im, while the size 
changes in layer 3 (the NH 4 SH layer) have the largest 
impact between 2.5 and 2.9 /rm. If we had used as the 
example a fit at low phase angles, and included an ap¬ 
propriately large number of quadrature points to capture 
the backward phase function peak, then the large parti¬ 
cle layer (layer 3) would also have an important effect 
at short wavelengths. But models with such phase func¬ 
tions don’t fit as well as models without strong backward 
peaks. 


In spite of the distinctive derivative spectra shown in 
Fig. [23l there are many spectra where the best fit NH 3 
and NH 4 SH layers are so close together that they be¬ 
come much more highly correlated with respect to pres¬ 
sure and optical depth variations. To fit these spectra 
and derive reasonable uncertainty estimates, we revised 
the model to use a fixed pressure offset between layers 
2 and 3, so that the uncertainty for the layer 2 pressure 
then becomes the same as that for the layer 3 pressure, 
except that the two layers have correlated uncertainties. 
This tendency for both small NHa-coated particles and 
much larger NH 4 SH particles to reside at nearly the same 
pressure level suggests that they might actually be mem¬ 
bers of two different populations in the same layer. It is 
also possible that NH 3 forms as a coating on the larger 
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Fig. 20. — Pressure vs. optical depth parameters for 4-layer 
multiple-scattering models that provide best fits to VIMS spectra, 
as given in Table [4| The pair of triple-dot-dash curves are cu¬ 
mulative optical depth profi les derived from belt and zone opacity 
profiles of lirwin et al.l 1120011') 



Lotitude, ° 



Lotitude, ° 

Fig. 21.— Optical depth vs. latitude (top) and pressure vs lati¬ 
tude (bottom) for the 4-layer model with NH3 in layer 2 (plotted 
as asterisks) and NH4SH in layer 3 (plotted as open circles). Layer 
1 pressure are excluded from the bottom plot because they are not 
very meaningful. 

NH 4 SH particles, although that would not work well for 
the entire NH 3 contribution, because the small particle 


component is needed to raise the reflectivity at short 
wavelengths, and the coated large particles don’t seem 
capable of producing the 2.97-/rm feature. 

5.6. Sensitivity of spectra to absorber gas profiles. 

The effects of different ammonia and para hydrogen 
profiles are illustrated in Fig. [24l The difference be¬ 
tween NEB and EZ profiles is relatively small, and our 
best fits were obtained with the the EZ profile for all 
spectra for which we did fits with both profiles, sug¬ 
gesting that the NEB prohle we chose did not contain 
quite enough NH 3 at high altitudes. The probe pro¬ 
hle makes a huge difference and is not a best ht with 
even the darkest spectrum we observed (13xl4y). While 
we cannot rule it out in a local region, due to lack of 
spatial resolution, it is certainly not acceptable over an 
extended region comparable to the area sampled by a 
VIMS pixel. The difference between normal and equi¬ 
librium hydrogen prohles is also relatively small and re¬ 
stricted to a region that does not provide a signihcant 
constraint on the cloud models. Given that atmosphere 
appears t o be intermediate between t hese states at low 
latitudes (jConrath and Gieraschlll98ll , the actual differ¬ 
ence is about half that between the extreme cases shown 
in the hgure. Thus, we did not try to account for po¬ 
tential variations in the ortho-para ratio. Instead we 
assumed equilibrium hydrogen in all the hts, conhdent 
that this convenience would have no significant effect on 
our conclusions. 

5.7. Effects of calibration scale-factor errors 

We fit several spectra with the nominal calibration and 
also with the VIMS I/F spectra multiplied by a scal¬ 
ing factor to simulate the effect of a change in abso¬ 
lute calibration. For a uniform fractional decrease in the 
measured I/F, we found that derived pressures and opti¬ 
cal depths also decreased, the former by relatively little, 
just a fraction of the fractional calibration change, and 
the latter by much more, changing by one (at low op¬ 
tical depths) to several times (at higher optical depths) 
the fractional change in calibration (large optical depths 
require a larger fractional change to produce the same 
fractional change in reflectivity). This is not too sur¬ 
prising, since the overall brightness (which is changed by 
the scaling) is a strong function of the aerosol loading, 
while the shape of the spectrum (which is not changed) 
is strongly controlled by the vertical distribution of the 
aerosols. These results suggest that opacities we derive 
have additional uncertainties contributed by calibration 
uncertainties of about 5% to 20%, while the pressures 
have additional uncertainties 0-4%. We also found about 
a 10 % affect on the wavelength dependence exponent for 
the bottom layer. The particle radius error contribution 
was not significant for layer 2 but 5-10% for layer 3. All 
these contributions should be combined with prior esti¬ 
mates using the root-sum-square method. 

6. DISCUSSION 

6 .1. Comparison with Galileo Probe results 

The Galileo Probe entered the atmosphere of Jupiter 
on 7 December 1995 at a latitude of 6 . 6 °N, near the 
boundary between the Equatorial Zone and NEB, and 
at the edge of a 5-^m hot spot, which is a region of un¬ 
usually low aerosol opacity (in the deeper layers). The 
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Fig. 22.— Best-fit spectrum for location llxl4y (thick gray curve) compared to spectra in which one layer is removed from the model, 
for each of the layers. Layer 2 is the cloud of NH 3 -coated particles and layer 3 is the NH 4 SH cloud. Note that no values were computed 
between 1.58 fim and 1.68 /rm. 


Nephelom eter (iR agent et al.||199^ and the Net Flux Ra¬ 
diometer (iSromovskv et al.lll998D both detected a cloud 
near the 1.2-bar level (base near 1.35 bars) at the edge of 
a hot spot, with roughly unit opacity at a wavelength of 
5 /tm. This pressure and optical depth are roughly com¬ 
parable to that of the bottom cloud derived from several 
of the VIMS spectra, but there is not a good match to 
any obtained close to the probe entry latitude. This mis¬ 
match is not too surprising; given the large size of the 
VIMS pixels and the small scale of a typical hot spot, a 
close agreement cannot be expected. 

An additional cloud layer was detected indirectly. A 
cloud above the level at which probe measurements be¬ 
gan (440 mb) was inferred from Net Flux Radiometer ob¬ 
servations of spin-induced modulations of the direct solar 
beam signal at multiple wavelengths (|Sromovskv et ahl 
I1998H . The optical depth was estimated to be ^1.5-2 at 
a wavelength of 0.5 /rm, and the particle size between 
0.5 and 0.75 /rm. This upper cloud seems reasonably 
consistent with our VIMS results for nearby latitudes 
( 12 xl 6 y and 26x31y). As clouds close to this level are 
found throughout the equatorial zone, it is plausible that 
these might also be present at the edge or even over a hot 
spot. Although we find a combination of small and large 
particles, their spectral effects on the solar beam inten¬ 
sity in the NFR passbands might be similar to those of 
a single population of particles of intermediate size. Our 
lowest opacity model (for the 12 xl 6 y spectrum) has two 
contributions that would come from above the probe de¬ 
ployment level. At 2.0 /rm the NH 4 SH layer has a small 
optical depth of 0.23, which would not increase much at 
0.5 /im because of the large particle size in this layer. 
The layer of NHa-coated particles, with an optical depth 
of 0.37 at 2 /im, would provide most of the opacity at 
0.5 /im, which could plausibly 4-5 times greater than its 


2-/im value. This would yield a total 0.5-/tm opacity 
comparable to the probe estimate, which is for a layer 
that does not vary much in optical depth for any of our 
spectral models, making it more likely that the probe 
comparison is a meaningful one. 

Our VIMS model results also agree with a lack of cloud 
particles in the 700 mb region, and we do find lower cloud 
layers close to the probe cloud density maximum near 1.2 
bars, although the optical depth we find there is generally 
much higher than obtained by the probe, which was in 
an unusually transparent region of the atmosphere. We 
don’t agree with the probe-derived ammonia profile how¬ 
ever; when we try to fit the least cloudy region (I3xI4y) 
with that ammonia profile, the best fit value is nearly 
double what we obtained using either equatorial zone 
or NEB ammonia profiles. The part of the profile that 
presents problems is the low level of NH 3 in the 400 mb 
to 1.5 bar pressure range, which is the part derived from 
NFR net flux measurements. 

6 . 2 . Comparison with Galileo NIMS results 

In Fig. [25] we compare our opacity results at a refer ¬ 
ence wavelength of 2 ttm to th ose of lirwin et al.l ()2nni[ ) 
and lirwin and Dvudii^ ()2002D . These comparisons are 
made in terms of the extinction optical depth per unit 
In P interval. This is roughly equivalent to optical depth 
per unit altitude. We estimated that function for our re¬ 
sults by convolving our opacities at discrete layers with a 
Gaussian smoothing function of the form exp {—x^/uP'), 
where the width w was set to roughly simulate the ver¬ 
tical cor relation used in th e NIMS retrievals. We didn’t 
scale the lirwin et al.l (| 200 lD results because they assumed 
wavelength-independent para meters over the entir e 1-2.5 
/im interval. We did scale the lirwin and Dvudinal (120021) 
results from their 1 -fim reference wavelength to our 2 -fiui 
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Fig. 23.— Best-fit model spectrum for location 25x35y (upper left), and fractional derivatives of that spectrum with respect to model 
parameters for pressure (PI, P2, P3, P4), optical depth (TAUl, TAU2, TAU3, TAU4), particle radius (R2, R3), and A exponent for layer 
4 (WLEXP4), in order from left to right and top to bottom. Here model layer 2 is made of synthetic NH 3 -coated particles and layer 3 of 
NH 4 SH. 


reference using their derived particle radius profile and 
calculations of the wavelength dependence of Mie scat¬ 
tering efficiency (using their refractive index of 1.4-|-0i). 
We find a rough agreement on the fractional increase 
of opacity with depth and agreement on the pressure 
level of the local peak in opacity, but don’t agree well 
on the absolute values of opacity. Our results contain 
several times as much optical depth as the obtained by 


lirwin et al.l (120011 ). and disagree e ven more with the pro- 
file of lirwin and Dvudinal (I2002D . which has nearly an 
order of magnitude less opacity when scaled to a wave- 
length of 2 u m. A t that wavelength it also disagrees with 
lirwin et al] (|2001l) . 

A factor that might i ntroduce d is agree ments 
with the NIMS analysi ^s of lirwin et al.l (I2001D and 
lirwin and Dvudinal (120021) is the difference in treatments 
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Fig. 24.— Spectra (top) for different gas opacity profiles using the same aerosol structure (the best fit structure for 25x35y and the EZ 
profile), and ratio (bottom) of each spectrum to the EZ spectrum. Note that no values were computed between 1.58 fim and 1.68 /im. 




Fig. 25.— Exti nction optical depth per unit per unit InP interval for a belt sample (left) and for a zone sample (right) derived by 
lirwin et al.l (120011) fsolid lines). The over-plotted circles are optical depths from VIMS fits for the 13xl4y spectrum (left) and the llxl4y 
spectrum (right), which represent belt and zone cloud structures respectively. The dashed lines are the dr/dlnP values derived from the 
smo othing of VIMS scal e d laye r optical depths, computed as described in the text. The dotted curves (left panel) is the opacity profile 
from ITrwin and Dvudin^ (120021) scaled to a wavelength of 2 fim, as described in the text. 


of NH 3 absorption. ISromovskv and FrvI (|2010f) show 
that prior ammonia absorption models provided too 
little absorption at 1.59, 1.9, and 2.3 fim. As can 
be seen from Fig. [211 variations in the amount of 
ammonia absorption have a significant impact on the 


I/F observed in the deeply penetrating windows at 1.28 
and 1.59 ^m, which provide important constraints on 
the vertical location of aerosols. Likewise, variations 
in NH 3 absorption models might affect the pressure 
derived from fitting the spectra. 
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Another difference that might affect the comparisons 
is that our fits result in different particle sizes, with the 
NH4SH lay er having much larger pa rticles than were 
derived by lirwin and Dvudiii^ l2002fl . Our NH3 and 
NH4SH layers ar e made with refracti v e inde xes of specific 
materials, while lirwin and Dvudin^ (j2002fl used conser¬ 
vative particles, have a strong influence on scaling results 
to other wavelengths. 

6 . 3 . Comparison with Galileo SSI results 

Our results are also in partial agreement with results 
based on Galil e o Solid State Imager (SSI) observations. 
iBanfield et al.l (Il998bll . who analyzed SSI data in three 
channels at 756 nm, 727 nm, and 889 nm, conclude that 
below an ubiquitous upper tropospheric haze, which ex¬ 
hibits little small scale variation, there is a highly variable 
cloud component (0 to 20 in optical depth), usually at p 
= 750±200 mb, less than a scale height in vertical extent, 
and the principal cause of features seen at red and longer 
wavelengths. At first, this appears to be in serious con¬ 
flict with our VIMS results, but it is only the lower end 
of the pressure range that is most in conflict, as a com¬ 
parison of cumulative opaci ties will soon demonstrate . 
The subsequent analysis by ISimon-Miller et ahl (|2001f) , 
reached similar conclusions, though somewhat less dis¬ 
crepant with near-IR results, specifically that a white 
cloud exists at pressures of 880-970 mb in belts and 620- 
720 mb in zones. These models use a relatively simple 
cloud structure with a small number of parameters, be¬ 
fitting the small number of spectral constraints provided 
by the observations. The structure consists of a strato¬ 
spheric haze extending to a pressure that marks the top 
of the main cloud layer with a bottom pressure at which 
(generally) is located a physically thin sheet cloud. The 
main cloud is assumed to have the same scale height as 
the gas, leaving two pressures and three optical depths 
to be constrained by the models. 

In Fig. [21] we provide comparisons for a belt (panel 
A) and a zone (panel B) betwe en our VIMS c umula- 
tive optical depths and those of IBanfield et al.l (|1998all 
and ISimon-Miller et al.l (120011) . Here our VIMS-derived 
opacities are scaled to their wavelength of 0.75 ^m. We 
scaled up the stratospheric haze opacity by a factor of 19 
(assuming a dependence) and the layer near 300 mb 
by a factor of ten, which is appropriate for our compos¬ 
ite particle with a radius of 0.3 pern. The bottom layer 
was scaled up by a factor of 2.67, to account for the 
approximate dependence we inferred for that layer. 
The main NH4SH 500-600 mb layer, was not scaled be¬ 
cause its opacity is relatively flat for the best-fit particle 
sizes, although its reflectivity can increase significantly 
at shorter wavelengths due to its high real index and re¬ 
duced absorption. At SSI wavelengths this layer appears 
to make a relatively small contribution compared to the 
layers surrounding it. 

Turning to the SSTVIMS comparison itself, it is not 
surprising that our stratospheric haze opacity is much 
too high because of the excessive I/F values in the VIMS 
spectra near 2.3 /rm . Otherwise, there is gross agreement 
on the opacity levels deeper in the atmosphere. For the 
belt comparison (panel A) the cumu lative opacity max¬ 
imum is in close agreement with the ISimon-Miller et al.l 
(|200If) result, while for the z one comparison (panel B) 
the best agreement is with the lBanfield et al.l (|1998bll re¬ 


sult, although there is a significant disagreement in the 
vertical location of the high-opacity bottom cloud. Our 
VIMS analysis puts it close to l.I bars, while the SSI 
results place it closer to 700 mb. On the other hand, 
we do find a significant cloud layer close to 700 mb in 
this region, and perhaps with a different scaling from 
2 /rm to 0.75 /rm the disagreement might be substan¬ 
tially reduced, suggesting that the particle size is smaller 
or the inferred wavelength dependence is greater than 
what we obtain from our best fit over the 1-3.2 /rm re¬ 
gion. Another point worth remembering is that we are 
not comparing a large statistical sample, and some of 
the differences observed might be attributable to spatial 
and temporal variations. In the belt region, for example 
we obtain somewhat better agreement with the 25x35y 
spectrum. Clearly, a larger sampling of variations would 
allow for a more meaningful comparison. 

Features of the SSI analysis that might contribute 
the disagreement in bottom cloud locations include 
changes made in t he methane absorption coefficients of 
iKarkosch^ (|1998ll , failing to include hydrogen collision- 
induced opacity from the third overtone band, which 
provides significant opacity in the 756-nm window, and 
uncertainties in the actu al transmissions of the G alileo 
SSI filters. According to ISimon-Miller et al.l (|2001ll . the 
756-nm absorption coefficient was increased by a fac¬ 
tor of 10 and the 727-nm coefficient from 0.17 to 0.22. 
The latest constraints on these coefficients, given by 
IKarkoschka and Tomaskol (|2009f) . have resulted in only 
a factor of two increase in the 756-nm coefficient and 
no change in the 727-nm coefficient. The two Galileo 
SSI analyses thus place more absorption at higher al¬ 
titudes than current best estimates would imply, and 
thus the pressures of the deeper cloud layers might 
be underestimated somewhat. Particle size differences 
are also a potential so urce of disagreements on opacity. 
IBanfield et al.l (Il998bll use a fixe d particle size of 0.2 /xm, 
while ISimon- Miller et al.l (|2001ll made particle size a fit¬ 
ted parameter, but didn’t report specific values because 
they were poorly constrained. 

6 . 4 . Comparison with \Matcheva et al.\ \200A) CIRS 

results 

The iMatcheva et al.l (120051 ) analysis of 2000 Gassini 
GIRS observations of Jupiter is unusual in relying on 
emission rather than reflected sunlight and in not re¬ 
quiring the use of band models, but instead relying on 
line-by-line calculations. This work also provides an im¬ 
portant constraint on cloud composition because they 
are sensitive primarily to absorption, whi ch varies greatly 
from o ne material to another (Fig. 1131) . IMatcheva et al.l 
([200511 chose a narrow spectral window centered at 7.18 
/im, which is free of NH3 absorption and thus indepen¬ 
dent of the NH3 mixing ratio profile. Their analysis ig¬ 
nores scattering, based on the assumption that all par¬ 
ticle radii are less than 1 /xm, although this is contra¬ 
dicted by other modeling effor ts, e. g. this work as we ll 
as that of IBrooke et al.l ( 19981) and IWong et al.l (1200411 . 
which generally find a significant component of large par¬ 
ticles (~10 /xm in radius) within the region relevant to 
the Matcheva et al. analysis. At latitudes with relatively 
low cloud opacity, such as the NEB, they find a cloud ab¬ 
sorption peak at 1100±50 mb, and a cumulative optical 
depth near 1 down to 1400 mb. At latitudes with thick 
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Fig. 2 6.— Cumulative optical depth for a belt sample (A ) and for a zone sample (B) derived from Galileo SSI observations bv IBantieid et all 
Ill998bl ) (dashed curves) and bv ISimon-Miller et al.l 1120011 1 (dot-dash curves). The over-plotted circles are optical depths from VIMS fits for 
the 13xl4y spectrum (A) and the llxl4y spectrum (B), which represent belt and zone cloud structures respectively. VIMS layer optical 
depths (circles) and cumulative optical depths (solid line s) have been scaled to a w avelength of 0.75 fim as described in the text. The 
triple-dot-dash curve is the cumulative opacity profile from lirwin and Dvudin^ II2002I) scaled to a wavelength of 0.75 fim, using their particle 
radius profile. 


overcast such as the EZ, they find peak absorptions at 
p> 900 mb and optical depths up to t > 4. Their tech¬ 
nique is limited to total optical depth (cloud plus gas) of 
r < 4. In the EZ, this limits retrievals to p < 980 mb, 
and in the NEB to p < 1.5 bars. 

In Fig. [27] we compare our absorption optical depth 
profiles inferr ed from the 13xl 4 y and llxl4y VIMS spec¬ 
tra with the iMatcheva et al.l ((20051) absorption optical 
depth profiles for a region near 12° N (left panel) and 
the equator (right panel), respectively. The comparison 
is made in terms of the parameter dTa/dhiP, which can 
be thought of as a crude measure of optical depth per 
unit height. Our model opacities from VIMS are scaled 
to a wavelength of 7.18 /im to make a more appropriate 
comparison to the Matcheva et al. results. The VIMS 
stratospheric opacities, scaled using a dependence, 
are too small to be relevant at 7.18 ^m. For the 300 mb 
layer we scaled from 2-/im values to 7.18-^m values us¬ 
ing the factor of 57 reduction that applies for our 0.3-^m 
composite particle, and a further reduction to account 
for the absorption optical depth being 80% of the ex¬ 
tinction value. These scalings make this layer ignorable 
as well. The optical depth of the main 3-/im absorbing 
layer has a relatively fiat wavelength dependence and was 
used without adjustment, except that we accounted for 
absorption optical depth being only 43% of the total. At 
7.18 /rm the particles in this layer have a single scatter¬ 
ing albedo of ^0.57 if they are pure NH 4 SH compared to 
~0.7-0.85 if they are pure NH 3 . 

The appropriate scaling for the VIMS bottom cloud is 
unclear. The wavelength exponents for the bottom layer 
in the reflecting layer fits (Table [T]) are generally small 
and slightly negative for most models. However, the ex¬ 
ponents obtained from the multiple scattering fits (Table 


a vary widely, being mostly negative for the low phase 
angle fits and mostly positive for the medium phase an¬ 
gle fits. We decided to use n=-I for the low opacity 
example and n =0 for the high-opacity equatorial exam¬ 
ple, mainly because they provided a crude agreement 
with the Matcheva et al. results at higher pressures. 
For this bottom cloud we also assumed the same absorp¬ 
tion to extinction optical depth ratio as for NH 4 SH. Af¬ 
ter scaling our opacities to the appropriate wavelength, 
we distributed them vertically using the same Gaussian 
function described for the NIMS comparison, this time 
adjusting the vertical width to simulate the vertical res¬ 
olution of the Matcheva et al. retrieval. This was done 
to facilitate comparison and not because we believe that 
the opacity is so distributed. 

At 12° N, the VIMS and CIRS opacity profiles are in 
crude agreement on vertical gradient, total absorption, 
pressure location of the absorption peak (near I.l bars), 
but disagree on the magnitude of the peak. This is prob¬ 
ably about as good an agreement as might be expected, 
given the very different models for vertical distribution. 
For th e 12° N case, we a l so inc luded a scaled version 
of the lirwin and Dvudinal (|2002l) profile, but added an 
imaginary index of 0.02 to simulate NH 3 absorption and 
0.2 to simulate NH 4 SH absorption. Here we see that the 
NIMS result does not scale well to the 7-/rm region, no 
matter which absorber we assume. The vertical gradient, 
peak location, and magnitude are in strong disagreement 
with the other results. 

The equatorial comparison in Fig. [23 tells a somewhat 
different story. Our results provide a weaker vertical 
opacity (after smoothing), while the CIRS analysis yields 
a very strong gradient. It is worth noting that the lay¬ 
ers in disagreement are those in which we find significant 
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scattering contributions, suggesting that the Matcheva 
et al. analysis might need to be revised to take account 
of scattering effects. The comparison might also be af¬ 
fected by local variations, and might be better if we had 
a larger statistical sample. It may also be a problem 
for one or the other analyses to obtain accurate results 
for regions with very high cloud opacity, since sensitiv¬ 
ity beyond the first few optical depths is very diminished 
for both VIMS and CIRS analyses. Perhaps the great¬ 
est uncertainty in this comparison is associated with the 
uncertainty in how to scale the bottom cloud results to 
the same wavelength. Currently we can only say that the 
two results might be consistent. 

6.5. Cloud composition 

Models that fit the VIMS spectra have some ammonia 
ice in the 330-550 mb region, probably in the form of 
ammonia combined with conservative material (near 3 
/rm ) of relatively high real refractive index. NH 3 could 
be present either as a coating or as the core material, in 
either case representing less than half the total particle 
volume. These particles have sub-micron particle radii, 
probably near 0.3 /im. The non-NHa component might 
be sedimented hydrocarbons generated by photolysis of 
stratospheric methane, although it does not appear likely 
that hydrazine is a major component. The main layer 
responsible for the 3-pia absorption feature is generally 
in the 500-600 mb region, and the best-fit composition 
of this layer is NH 4 SH, although this layer could also 
contain a contribution from NH 3 ice, either in the form 
of a separate particle population or perhaps as a coating. 
In fact, the close proximity of the NH 3 and NH 4 SH layers 
in many cases, suggest that they might be co-mingled. 

It is possible within the constraints of the spectra that 
low-opacity regions have 3-^m absorption provided en¬ 
tirely by NH 3 and that NH 4 SH is only needed in the 
higher-opacity regions. But it seems more sensible that 
the deeper clouds are all composed of similar materials, 
and that NH 3 contributions appear in a different layer 
at generally lower pressures, although the fact that the 
two layers are virtually on top of each other in the low 
opacity regions raises questions about that possibility as 
well. Perhaps the entire particle distribution is coated 
with NH 3 , but that only the small-particle component 
actually allows the NH 3 absorption features to become 
visible. Perhaps the top layer of clouds is the small- 
particle component that displays these features, while 
the bulk of the cloud is larger particles that don’t show 
the features. 

We tried a number of models for the bottom cloud. The 
one that seemed to work best is one in which the phase 
function is wavelength ind ependent and matches that of 
ISromovskv and Fr^(|2002ll . while the optical depth varies 
as A", with n ranging from -1.3 to -I-1. 6 . This resulted 
in pressures from 800 mb to 1.3 bars and optical depths 
from 5 to 20 (Table S]). This cloud might be composed 
of NH 4 SH, but using that material in spherical parti¬ 
cles degraded the fit quality. If we used a large number 
of quadrature points to capture the wavelength depen¬ 
dence of the phase function with higher accuracy, we 
found the details of that variation did not match the 
observations. On the other hand, if we used a small 
number of quadrature points to suppress the backscatter 
peak and its variations, this also made the lower cloud 


insufficiently bright, resulting in poor fits to the window 
regions. It appears that if NH 4 SH is the main compo¬ 
nent of this cloud, it must have a lower imaginary in¬ 
dex in the near-IR window regions, or have increased 
backscatter with a wavelength dependence different from 
that of spherical particles. There may also be accompa¬ 
nying variations in deeper regions of the atmosphere that 
are not constrained by the observations we used. 

6 .6. Dynamical considerations 

At first glance, a scenario in which NH 4 SH provides the 
main component of cloud particles in vicinity of 500 mb 
seems implausible. Such a scenario is certainly not con¬ 
sistent with a uniform upwelling in zones and descend¬ 
ing air in belts, as common ly suggested in early stud ¬ 
ies of Jovian dynamics (e.g. I Hess and Panofsk"^ (II969I1 : 
IStonel (IIOTGIH . However, there is an observational con¬ 
flict with that simple belt-zone dynamical model, namely 
the Galileo discovery that most of the lightning occurs 
within belts, where air was thought to be descending, and 
is associated with rapidly growing cloud clusters reach¬ 
ing p ressures of a few hundred millibars (|Ingersoll et al.l 
120001) . These events would certainly dredge up conden¬ 
sible material from deeper levels, including NH 3 , H 2 S, 
water, and likely particles of NH 4 SH. If these events 
provided a major source of condensibles in the upper 
troposphere, then it is less difficult to understand how 
NH 4 SH particles might contribute a major source of 
cloud material near 500 mb. The potentially important 
role of these thunderstorms in controlling the compo- 
sition of the upper troposphe re is also highlighted by 
[Showman and de Pateij (120051) . who provide a dynam¬ 
ical explanation for another non-intuitive observation, 
namely the global depletion of NH 3 vapor between 500 
mb and 5 bars in both belts and zones. Their explana¬ 
tion is based on the idea “( 1 ) that the majority of air 
that ascends across the 5 bar interface resides in isolated 
thunderstorms, and (2) that 50% of the ammonia within 
these storms is lost through the 4-6 bar level either by 
direct ammonia rainout through the base of the storms 
or downward transport of ammonia vapor in convective 
downdrafts moistened by evaporation of rainfall.” They 
speculate that “air undergoing large scale ascent above 
the 0.5-bar cloud tops in zones is probably supplied hor¬ 
izontally from the belts (where thunderstorms predomi¬ 
nantly occur) rather than from below.” T he belt-to-zone 
mass flux could be provided by eddies (|Ingersoll et al.l 
l2000() . It also seems plausible that during the strong con¬ 
vective events, composite particles could form, including 
ammonia-coated NH 4 SH. When distributed horizontally, 
such particles would be mixed with less ammonia-rich 
air, promoting evaporation of the coating and leaving 
behind the NH 4 SH core, which might also contain some 
component of H 2 O, although the VIMS spectra would 
not be consistent with a large fraction of water ice in 
such composite particles. A quantitative microphysical 
investigation of this scenario is needed to provide a better 
understanding of its plausibility. 

6.7. Relation to ammonia ice detected at thermal 
wavelengths 

Detections of N H 3 spectral f e ature s near 10 /rm were 
first obtained by iWong et al.l (120041) . They used the 
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Fig. 27.— Cloud absorption optical depth per unit InP interval for 12° N (left) and the equator (right) derived bv IMatcheva et'a n 11200511 
from CIRS thermal emission observations at 7.18 fim. Above and below the plotted pressures, the retrieval is dominated by guess profiles 
and thus not shown. The over-plotted circles are optical depths from VIMS fits for the 13xl4y spectrum (left) and the llxl4y spectrum 
(right), which represent belt and zone cloud structures respectively. The scaling of VIMS results to 7.18 fim is described in the text. The 
dashed lines are the VIMS results con volved with a smoothing f unction that crudely simulates the vertical resolution of the Matcheval et 
al. retrieval. The dTa/d\nP values for lirwin and Dvudinal II2002I 1 (dotted lines at left) are scaled according to particle size using imaginary 
indexes of NH 3 {Ni=0.02) and NH 4 SH (Afi=0.2) to simulate NH 3 and NH 4 SH absorption respectively. 


brightness temperature difference between 1040 cm“^ 
and 1060 cm“^ to measure the strength of the ammo¬ 
nia ice signature. This parameter peaked at the equator 
and at 23° N, and was much more widely distributed 
than t he very rare discret e ammonia ice clouds identi¬ 
fied bv I Baines et al.l (| 2 nn 2 f l. but was clearly not as om¬ 
nipresent as the 3-pm absorption feature. They modeled 
their observations with 4:1 prolate spheroids of a volume 
equivalent radius of 0.79 pm, a cloud base of 790 mb, 
and extending up to 100 mb with a 1:1 particle to gas 
scale height, accompanied by a much more compact ( 1/8 
the gas scale height) gray absorber, which they modeled 
as 10-pm NH 3 ice spheres. For average spectra from two 
regions (22-25° N, covering 140-240° W, and 14-17° N, 
covering 10-70° W) they obtained for the small parti¬ 
cle component optical depths of 0.75 and an upper limit 
of 0.2 respectively. The larger optical depth is likely 
roughly comparable to optical depths in the equatorial 
region where they also found a peak in the ice signature. 

Although the v ertical variation in opacity of the 
iWong et al.l (|2004D cloud models is not a good match 
to ours, the latitudinal variation in ammonia ice they in¬ 
fer is at least qualitatively consistent with our results in 
Fig. EH In our models the NHa-containing layer descends 
from 400 mb to 600 mb from the equator t o 14°N. This 
chang e is especially relevant because of the IWong et al.l 
(|2004l) conclusion that ammonia ice present deeper than 
500 mb would not be detectable at 10 pm because it 
would be too far below the peak of the contribution func¬ 
tion of the observations. This implies that their obser¬ 
vations would detect much less NH 3 in the NEB. Their 
broad vertical distribution of the small particle compo¬ 
nent for the cloudiest region, which should be similar to 
the equatorial region, would place about half its opti¬ 


cal depth at pressures less than 450 mb, which would 
be about 0.38 extinction optical depths, comparable to 
what we would estimate for the layer of small ammonia- 
coated particles. Their more compact cloud could easily 
be composed of NH 4 SH. 

6 . 8 . Relation to SIACs 

iBaines et al.l (| 2002 D identified anomalous discrete 
cloud features in the wake of the GRS, which they termed 
spectrally identifiable ammonia cloud (SIACs) based on 
low reflectivity at 2.7 pm relative to high reflectivity at 
1.6 pm, and a local dip in reflectivity at 2 pm relative 
to 1.94 pm, both depressions characteristic of ammonia 
ice absorption (the longer wavelength absorption is only 
apparent in optically thick clouds, and not noticeable in 
Fig. fTH) . These features were observed in Galileo NIMS 
spectra, taken at very high spatial resolution, and appear 
to occupy less than 1% of Jupiter’s cloud-top surface. 
The VIMS spectra we analyzed (or any other VIMS spec¬ 
tra) have far too low a spatial resolution to detect SIACs. 
Unlike SIACs, the 3-pm absorption feature we analyzed 
appears to be very widely distributed on Jupiter. Al¬ 
though we focused on low-latitude spectra for detailed 
analysis in this paper, we did look at VIMS spectra over 
a wide range of latitudes and longitudes, and found ob¬ 
vious strong absorption at 3 pm everywhere we looked. 
The SIACs appear to be the result of unusually strong lo¬ 
calized vertical transport that causes significant conden¬ 
sation of ammonia vapor, such that the ammonia spectral 
features become significantly enhanced. It appears from 
our results, that most locations on Jupiter don’t have 
very much condensed NH 3 , so that the the spectral char¬ 
acter of SIACs is rarely observed. The alternative possi¬ 
bility that NH 3 is in fact the major cloud component, but 
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has its spectral signature significantly altered by some 
photochemical “tanning” process remains to be quanti¬ 
tatively investigated, and whether it could satisfy known 
spectral and photochemical constraints seems doubtful. 

6.9. Comparison with Saturn 

Saturn’s near-IR spectrum does not contain the broad 
and strong 3-/im absorption feature so apparent in 
Jupit er’s spectrum (lEncrenaz et al.l 119991 : iBaines et al.l 
l2005f) . The stark difference between spectra of these 
two planets is illustrated in Fig. 1281 where we compare a 
typical low-latitude Saturn spectrum from VIMS data 
cube CM_1587635989 with a VIMS Jupiter spectrum 
(12xl6y). The high reflectance of Saturn’s clouds at 3.0 
^m (~0.2), just where Jupiter’s I/F is quite low (^0.02), 
implies that the main cloud layer contributing to Sat¬ 
urn’s reflection spectrum in the near IR is not composed 
of medium sized pure ammonia ice particles (nor NH 4 SH 
particles) . This is likely the layer of upper troposheric 
haze that iPerez-Hovos et al.l (|2005[ 1 found to extend from 
a base at 300 mb - 500 mb to a top near 40 mb - 80 mb, 
varying with latitude and season, but generally of sub¬ 
stantial optical depth 6-20 at ^1 /rm), with a weak 
dependence on wavelength. This layer is probably thick 
enough at 3 /rm to obscure the much deeper ammonia 
cloud thought to reside between 1.4 and 1.8 bars. 



Fig. 28.— A typical low-latitude VIMS spectrum of Saturn (from 
x=21, y=2 in data cube CM_1587635989) reveals absorption by 
PH 3 between 2.8 pm and 3.0 pm, but does not display the broad 
absorption from 3.0 pm to 3.2 pm that is characteristic of the 
typical Jovian spectrum, such as the plotted 12xl6y VIMS spec¬ 
trum (dashed). The cloud albedo (dot-dash) inferred bv IKim et al.l 
ll 200 g'l from analysis of a groundbased high-resolution spectrum 
contains a local absorption qualitatively consistent with a layer of 
small particles containing a minor fraction of ammonia ice, either 
as a core or coating. The vertical dotted line is at 2.965 pm. 

However, there is some evidence for at least a compo¬ 
nent of ammonia ice in Saturn’s tropospheric haze. From 
analysi s of a high - resolu tion groundbased spectrum of 
Saturn, iKim et ahl (120061 ) developed evidence for a weak 
cloud absorption feature at 2.96 /tm, where their inferred 
model cloud layer at 460 mb had its albedo drop by 
^ 20 %, which they suggested might be due to ammonia 
ice. A small absorption feature at the same wavelength 
is also seen in many spatially resolved VIMS spectra of 
Saturn (also evident in Fig. [28l), and is not a character¬ 
istic of the phosphine spectrum that provides the main 
absorption in this region. This very well may be due to a 
population of ammonia-coated particles within the main 
tropospheric haze layer of Saturn. A small composite 


particle of this type with ~30% ammonia content, would 
lack the broad absorption feature of larger pure ammonia 
ice, and have a relatively narrow, but subdued, peak ab¬ 
sorption displaced towards the ammonia imaginary index 
peak at 2.965 /rm, as is shown in Fig.[T3 This shifted and 
attenuated feature could also be produced by particles 
with a small core of ammonia ice coated by a thick layer 
of hydroca rbons, as demonstrated b y the benzene-coated 
example of lKaloeerakis et all (|2008D . It seems likely, but 
remains to be demonstrated, that a moderately thick tro¬ 
pospheric haze containing a component of such particles 
could provide accurate fits to Saturn’s near-IR spectra. 

7. SUMMARY AND CONCLUSIONS 

From an analysis of VIMS low and medium phase 
angle spatially resolved spectra obtained during the 
Cassini 2000 flyby of Jupiter, we obtained the following 
main conclusions regarding the VIMS instrument per¬ 
formance, the source of Jupiter’s 3-/im absorption, and 
Jupiter’s low-latitude cloud structure. 

1. By comparing VIMS spectra with groundbased 
spectra, selected NICMOS observations, and an 
ISO spectrum, we were able to confirm that the 
VIMS absolute calibration is relatively consistent 
with other observations, but that VIMS measure¬ 
ments have two significant artifacts, one in the 
1.60-1.68 ^m region where a huge responsivity cor¬ 
rection turns out to be insufficiently accurate, and 
the other is an excessively high I/F measurements 
at low signal levels, which caused model fits to yield 
stratospheric haze opacity ~10 times that derived 
from groundbased spectra . 

2. We found that VIMS spectra contain evidence 
for a prominent and widespread 3-/im absorber in 
Jupiter’s clouds and a much weaker absorber active 
near 2 /im. This was made obvious by reflecting 
layer fits constrained outside the absorbing region, 
which produce more than the observed reflectivity 
in the absorbing regions, by ^20% at 2.04 pm and 
by 300-500% at 3-3.1 pm . 

3. We were able to fit a low-opacity VIMS spectrum 
quite well using NH 3 ice as the only 3-pm absorber, 
but were not able to fit very well the spectra from 
much cloudier regions, the main defect being an 
absorption feature near 2 ^m in the model spec¬ 
tra that was much larger than observed. Models 
with NH 4 SH as the only absorber in a similar layer 
structure were more successful in fitting high opac¬ 
ity spectra and almost as good fitting low opacity 
spectra, the main defect being lack of a 2.97-^m 
spectral feature which appears to require a contri¬ 
bution from NH 3 . 

4. The best hts to VIMS spectra were obtained 
with both NH 3 and NH 4 SH absorbers present in 
Jupiter’s clouds, but the NH 3 contribution appears 
to help most in the form of a coating or as a core 
material comprising less than half the total particle 
volume. We roughly simulated NHs-coated parti¬ 
cles using a synthetic refractive index with absorp¬ 
tion similar to that of NH 3 but significantly weaker. 
A layer of these small (r^O.3 pm) particles above 
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a layer of much larger (r^lO ^m) NH 4 SH particles 
provided the best fits to all the spectra. The NH 3 - 
coated cloud layer (modeled as physically thin but 
ranging from 330 to 550 mb) was found to have the 
least variability in optical depth among the eight 
low latitude spectra, with an average of 0.38 and 
a standard deviation of only 0.05. The main 3- 
/rm absorbing layer varies in optical depth from 
0.23 to 2.1, a factor of nine, with pressures rang¬ 
ing from 358 to 635 mb. In most cases the NH 3 - 
containing layer was found in such close proximity 
to the NH 4 SH layer that they could be intermin¬ 
gled. A caveat in these results is that the simu¬ 
lation of the NH 3 -coated particle does not closely 
follow the size and wavelength dependence of a true 
NH 3 -coated particle, and thus some other structure 
may be required to provide physically consistent 
fits. 

5. Comparison with Galileo Probe results provided 
crude agreement on the vertical location of clouds, 
the upper level cloud opacity, the pressure of the 
lower cloud, and the relative lack of cloud parti¬ 
cles near 700 mb, but did not agree with the NFR- 
derived vertical profile of NH 3 . 

6 . Comparison w ith the analysis of Galileo NIMS ob¬ 
servations by lirwin et ahl (l 2001 f ) shows a rough 
agreement with the vertical trend of opacity and 
the location of the peaks in opacity, but we find 
greater absolute opacity values. We cannot con¬ 
firm (or rule out) the large belt-zone differences 
they fo und in the 1.5-2 bar regio n. When the re¬ 
sults of lirwin and Dvudinal (120021 1. which apply to 
belt regions, are scaled to 2 gm, there i s an e ven 
larger discrepancy than with lirwin et al .1 (l 2001 fl . 

7. Comparison of our results scaled to a wavelength of 
0.75 u,m wit h Galileo SSI results bviBanfield et al.1 
(Il998br) and iSimon-Miller et ahl (|2001h finds rea- 
sonable agreement with cumulative opacity pro¬ 
files, with the exception of the location of the lower 
high-opacity cloud in zone models, which the SSI 
analysis places near 700 mb and our analysis places 
near 1100 mb. 

8 . Comparison of our absorption profiles with 
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the analysis o f Cas sini CIRS observations by 
iMatcheva et al.1 (I2005D , with plausible scalings with 
wavelength, reveals a rough agreement, but with 
CIRS opacities less than ours by a factor of two or 
more. However, the validity of the comparison is 
uncertain because of the large uncertainty in ap¬ 
propriate scaling factors to use for the lower cloud. 

9. The presence of a layer of sub-micron NH 3 ice 
coated particles at pressures of 300-500 mb does 
not appear to conflict with thermal spectra in the 
9-10 fim region, where only a weak ammonia sig¬ 
nature is observed. This is possible because the 
coated particle can produce the reflective contri¬ 
bution needed to match the near-IR spectra, but 
ten times less extinction optical depth at 9.4 /im 
than would be obtained from pure NH 3 particles 
of the same size. Our estimated extinction optical 
depth for this layer is cr udely compatible w ith the 
optical depth estimate of lWong et al.1 (j2004[) in the 
latitude band we analyzed. The increased pressure 
we find for the NH 3 layer in the NEB is qualita¬ 
tively consistent with the low signature found by 
Wong et al. in a similar region. 

10. Saturn’s spectra do not exhibit the broad 3-/im 
absorption feature characteristic of Jovian spectra, 
perhaps because the cloud layer containing that ab¬ 
sorber is obscured by the moderately thick tropo¬ 
spheric haze that overlies it. However, Saturn spec¬ 
tra do exhibit a small absorption feature near 2.965 
/im, which is qualitatively consistent with NH 3 ab¬ 
sorption when diluted by combining NH 3 with a 
substantial conservative core or shell. 
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